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I .  INTRODUCTION 


i  .  Map  Projections 

The  subject  of  map  projections  has  developed  from  its 
inception  sometime  early  in  the  history  of  mankind,  consisting  of 
perhaps  a  rough  sketch  drawn  in  the  sand,  to  a  specialized  branch  of 
applied  mathematics  which  it  is  today.  Due  both  to  this  long  develop¬ 
ment  and  to  the  extent  of  the  subject,  it  is  possible  to  encounter 
numerous  definitions  of  a  map  projection.  The  very  name  of  the 
subject,  mnp  projections,  is  deceiving,  for  most  map  projections  are 
not  projections  in  the  true  mathematical  sense  of  the  vcrd.  In  the 
true  mathematical  sense,  the  term  projection  hus  the  connotation 
perspective  projection.  Perspective  projection  may  be  simply  described 
as  the  determining  of  the  position  on  it  pl  'ne  of  points  in  spice  as 
they  would  appear  on  the  plane  from  some  particular  fixed  view  point. 
Applying  the  above  uefinition  to  map  projection  would  indicate  that  a 
map  projection  is  the  representation  of  the  surface  of  "-he  earth,  or 
an  iccepted  datum  surface  approximating  the  surface  of  the  earth,  on 
a  plane  as  it  would  appear  on  the  plane  from  some  particular  fixed 
view  point.  However,  the  above  definition  can  certainly  not  be  used 
to  totally  describe  map  projection  in  the  sense  in  which  the  term  is 
commonly  used.  The  term  projection,  as  applied  to  map  projection, 
includes  numerous  geometrical,  mathematic al, ana  semi -geometrical 
transformations  of  a  iatura  surface  onto  a  plane.  Only  a  very  few  of 
these  projections  are  truly  perspective  projections  in  n  mathematical 


sense . 
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As  was  stated  earlier,  a  map  projection  has  oeen  detinea  in  many  ways. 
One  method  of  describing  n  map  projection  is  to  describe  the  trans¬ 
formation  of  the  meridians  and  para  lie  ls  irotn  tne  surface  of  the  earth, 
or  a  datum  surface  approximating  the  surface  of  the  earth,  to  a  map¬ 
ping  plane.  The  author  will  define  a  man  projection  as  a  aystematic 
arrangement  of  intersecting  lines  on  a  pl^ne  that  represent,  and  have 
a  one-to-one  correspondence  to,  the  meridians  and  parallels  of  the 
datum  surface.  These  lines,  drawn  to  represent  the  meridians  and 
parallels,  are  drawn  according  to  some  consistent  principle  in  order 
to  fulfill  certain  required  conditions.  Each  set  of  new  conditions 
produces  a  different  map  projection,  and  hence  there  potertially 
exists  an  unlimited  number  of  map  projections .  Common  conditions  used 
in  map  projections  are  such  that  distances,  areas,  or  angles  measured 
on  the  datum  surface  will  be  equal  to  those  measured  on  the  mapping 
plane.  In  addition  to  the  more  common  conditions  stated  above,  a 
myriad  01  other  conditions  have  been  used  to  produce  map  projections. 
The  conditions  by  which  each  particular  mnp  projection  is  designed 
are  naturally  b°sed  on  the  intended  purpose  of  the  map. 

1.2  Perspective  Map  Projection 
Based  on  the  discussion  in  Section  1.1,  it  is  apparent  that 
the  perspective  map  projection  is  only  one  of  the  many  different 
types  of  map  projections.  It  has  characteristics,  both  desirable  and 
undesirable,  which  dictate  its  use  or  non-use  for  particular 
purposes.  A  perspective  map  projection  may  be  ideally  vizualized  as  a 
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photograph  taken  of  a  portion  of  the  datum  surface.  A  photograph 
taken  by  an  aircraft,  a  rocket,  or  a  space  vehicle  of  the  earth,  moon, 
or  other  heavenly  body  is  actually  a  type  of  perspective  projection. 
This  certainly  is  one  of  the  interesting  aspects  of  the  study  of 
perspective  projections. 

Perspective  projections  are  based  on  simple  geometric 
principles.  A  straight  line,  vhich  will  be  called  the  projection 
axis  of  the  projection,  extends  from  a  point  on  the  d-itum  surface, 
through  a  point  in  space  named  the  projection  center,  to  the 
projection  plane.  This  projection  axis  is  perpendicular  both  to  the 
datut  surface  and  the  projection  plane.  The  ini-prsection  of  the 
projection  axis  and  the  projection  plane  locates  *>  point  denoted  as 
the  image  center.  Likewise,  all  other  points  on  the  datura  surface 
arf  located  on  the  projection  plane  by  extending  straight  lines  from 
the  particular  point,  through  the  projection  center,  onto  the  pro¬ 
jection  plane.  In  Figure  1,  C  represents  the  projection  centn-r  and 
C',  the  image  center.  P  is  an  arbitrary  point  on  the  latum  surface 
and  P'  is  the  pro  lection  of  this  point  on  the  projection  plane. 
Principles  of  elementary  geometry  prove  that  a  parallel  displacement 
of  the  projection  plane  along  the  projection  axis  chances  only  the 
scale  of  the  projection.  However,  a  change  in  the  location  of  the 
projection  center  will  chance  the  form  of  the  projection.  The  location 
of  the  projection  center  is  the  factor  vhich  determines  the  character¬ 
istics  of  a  particular  perspective  projection. 
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Figure  1 

One  means  of  classifying  perspective  projections  is 
according  to  the  orientation  of  the  projection  axis.  If  the 
projection  axis  is  coincident  with  the  axis  of  rotation  or  the 
datum  surface,  the  projection  is  referred  to  as  a  normal  or  polar 
projection.  If  the  projection  axis  lies  in  the  equatorial  plane 
of  the  datum  iurface,  the  projection  is  referred  to  as  a  transverse 
or  equatorial  projection.  If  the  projection  falls  into  neither  of 
the  above  two  special  categories,  then  it  is  referred  to  as  an 
oblique  projection.  It  should  be  kept  in  mind  that  the  normal  and 
transverse  projections  ^re  merely  special  cases  of  the  oblique 
projection  and  equations  lerived  for  oblique  projections  are 
applicable  to  all  special  cases. 
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1.3  The  Investigation  Of  Perspective  Map  Projections 
The  purpose  of  this  thesis  is  to  investigate  the  perspective 
projection  in  a  very  basic  manner.  The  primary  goal  is  to  develop 
general  mapping  equations  for  the  perspective  projection  suitable 
both  for  a  spherical  and  an  ellipsoidal  datum  surface.  The  term, 
mapping  equations,  refers  to  equations  which  transform  the  co¬ 
ordinates  of  a  point  from  a  latitude-longitude-elevation  three 
dimensional  type  coordinate  system  of  the  datum  surface  to  a  X-Y  plane 
coordinate  system  of  the  projection  plane.  A  discussion  of  the 
distortions  Inherent  in  perspective  projections  is  Included  to 

i 

illustrate  a  method  of  using  distortions  to  describe  a  projection. 

It  will  be  seen  that  a  combination  coordinate-distortion  type  table 
or  grid  is  an  effective  aid  in  visualizing  a  projection. 

The  investigation  is  divided  into  five  primary  parts.  In 
order  to  provide  general  background  material  and  also  provide 
equations  on  which  to  base  the  validity  of  equations  derived  later, 
Chapter  2  will  briefly  discuss  and  derive  in  a  conventional  manner 
the  mapping  equations  for  the  three  most  common  types  of  perspective 
projections.  These  three  common  projections  are  the  Gnomonlc 
Projection,  the  Stereographic  Projection,  and  the  Orthographic 
Projection.  Chapter  3  discusses  «  descriptive  geometry  approach  to 

i 

perspective  map  projections  as  presented  by  Erwin  Schmid  [12]. 

It  is  from  this  publication  that  the  author  first  became  interested 
in  the  subject  matter  of  this  thesis.  Chapter  4  is  concerned  with 
the  actual  development  of  general  mapping  equations  for  the 
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perspective  projection.  These  equations  are  developed  utilizing  a 
photogrammetric  approach,  an  approach  which  to  the  author’s  know¬ 
ledge  is  novel  to  the  study  of  map  projections.  Equations  are  first 
derived  using  the  sphere  as  a  datum  surface  in  order  "both  that  the 
validity  of  the  equations  may  be  checked  by  comparison  with  those 
equatious  derived  earlier  and  to  facilitate  for  the  render  the 
understanding  of  later  derivations.  In  an  analogous  manner  equations 
are  then  derived  using  the  ellipsoid  as  a  datum  surface.  It  is 
shown  that  the  mapping  equations  derived  for  an  ellipsoidal  datum 
surface  are  truly  General  Perspective  Projection  Mapping  Equations, 
suitable  for  both  spherical  and  ellipsoidal  datum  surfaces. 

Chapter  5  briefly  discusses  an  empirical  method  to  determine 
the  distortions  inherent  in  the  perspective  projection  of  a  latitude 
longitude  grid  of  a  datum  surface  onto  a  projection  plane.  Chapter 
6  illustrates  examples  of  six  types  of  perspective  projections  by 
means  of  tables  or  grids  containing  both  X-Y  plane  coordinates 
and  distortions.  The  coordinates  and  distortions  are  calculated  for 
the  intersection  points  of  latitude  and  longitude  lines  spaced  ten 
degrees  apart.  It  is  believed  that  the  examples  included  in  Chapter 
6  both  aid  in  the  vizualization  of  perspective  projections  and  at 
least  partially  validate  the  equations  derived  in  Chapter  4.  The 
computer  program  used  to  produce  the  examples  of  Chapter  C 
and  a  brief  explanation  of  the  computer  program  are  included  in  the 
Appendix. 
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2.  A  CONVENTIONAL  APPROACH  TO  PERSPECTIVE  MAP  PROJECTIONS 


Three  of  the  more  common  perspective  map  projections  are 
the  Gnomonic  Projection,  the  Stereographic  Projection,  and  the 
Orthographic  Projection.  The  mapping  equations  for  these  three 
projections  are  brief ly  derived  in  this  chapter  for  two  main  reasons. 
First,  the  derivations  will  illustrate  what  the  author  chooses  to  call 
a  conventional  approach  to  the  derivation  of  perspective  mapping 
equations,  similar  to  the  approach  used  in  []*].  Second,  the  equations 
derived  in  this  chapter  will  he  a  L«8is  of  comparison  for  those 
equations  derived  later. 


2.1  The  Gnomonic  Projection 

The  Gnomonic  Projection  is  developed  by  placing  a  projection 
center  at  the  center  of  a  sphere  and  projecting  this  sphere  onto 
a  projection  plane  tangent  to  it.  Figure  2  represents  a  sphere  with 
a  radius  R  and  a  center  at  C,  which  in  the  Gnomonic  Projection  is 
also  the  projection  center.  The  projection  plane  AB  is  made  tangent 
to  the  sphere  at  O'  and  this  point  is  the  image  center  of  the  pro¬ 
jection,  O'  is  located  at  a  latitude  #con  the  sphere.  P  is  an 
arbitrary  point  located  on  the  sphere  at  a  latitude  4>  and  at  a 
longitude  difference  of  X  from  O’.  P'  is  the  projection  of  P  on  the 
projection  plane. 

Figure  2  illustrates  the  following; 

CO'  -  CP  -  R 
Angle  O'NP  ■  X 
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Angle  NO'P  -  90  -  B 
Arc  NO*  -  90  -  * 

c 

Arc  NP  -  90  -  • 

Figure  3  1*  drown  In  the  plane  of  the  triangle  CO'P  and 
illustrates  the  projection  of  the  arc  O'P  onto  the,  projection  plane 
AB.  From  Figure  3>  it  is  evident  that 

(1)  O’P'  -  R  tan  q 

To  reduce  equation  (l)  to  rectangular  coordinates  on  the 
projection  plane,  the  coordinate  system  illustrated  in  Figure  4  Is 
developed.  In  this  coordinate  system,  O'Y*  represents  the  projection 
of  the  meridian  through  O'  and  O'X'  represents  the  projection  of  the 
great  circle  through  O'  which  Is  perpendicular  to  the  meridian 
through  O'. 

From  Figure  4,  the  following  is  derived: 

(2)  X’  .  ?P*  -  O'P'  cos  B-  R  tan  q  cos  B»  -  8l" 

C  O©  Q 

_  R  sin  a  sin  B 

(3)  Y' -EP' -O'P'  sin  B-R  tan  q  sin  B- - - 

Utilizing  the  properties  of  the  spherical  triangle  in  Figure  1,  the 

following  may  be  derived  by  the  Lav  of  Sines: 

sin  A  _sin  (90  -  B)  m  cos  B 
sin  q  sin  (90  -  #}  cos  § 

and  therefore 

(4)  sin  q  cos  B-  sin  X  cos  # 

By  the  Lav  of  Cosines,  the  following  is  derived: 

(5)  sin  q  sin  B  -  cos  |  sin  #  -  sin  #  cos  f  cos  X 

c  c 

(6)  cos  q  asin  f  sin  t+cos  I  cos  f  cos  X 

c  c 
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By  substituting  equations  (k) ,  (5),  and  (6)  Into  equations 
(2)  and  (3),  the  mapping  equations  for  the  Gnomonlc  Projection  are 
obtained: 


(7) 

(8) 


X'- 


y 


_ R  cos  f  sin  X 

sin  #„  sin  4  ♦  cos  9  cos  $  cos  X 

R  (cos  f  sin  f  -  sin  |  cos  f  cos  X) 

—i  . .  I.  ^ i  in  — a .Cl.  .  ■ 

sin  t  sin  |  ♦  cos  |  cos  ff  cos  X 
c  c 


2,2  The  Stereographic  Projection 
The  Stereographic  Projection  Is  developed  by  placing  a 
projection  center  at  a  point  on  the  sphere  vhich  Is  diametrically 
opposite  the  point  at  vhich  the  ^pjectlon  plane  Is  tangent  to  the 
sphere.  The  mapping  equations  for  the  Stereographic  Projection  are 
developed  In  an  analogous  manner  to  the  development  of  the  mapping 
equations  for  the  Qnoaonlc  Projection.  Figure  5  represents  a  sphere 
vlth  a  radius  R  and  a  center  at  0.  The  projection  center  Is  located 
at  C  and  the  projection  plane  AB  Is  tangent  to  the  sphere  at  O' .  O' 

Is  also  the  Image  center  of  the  projection.  P  Is  an  arbitrary  point 
located  on  the  sphere  at  a  latitude  #  and  at  a  longitude  difference  of 
X  from  O'.  P'  Is  the  projection  of  P  on  the  projection  plane. 

Figure  3  Illustrates  the  follovlng: 


OC  -  00'  -  R 
Angle  O'HP-  X 
Angle  MO'P.  90-6 
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Arc  HO'  ■  90  -  # 

c 

Arc  HP  »  90  -  # 

Figure  6  1*  dravn  in  the  plane  of  the  triangle  CO'P  and 
Illustrate!  the  projection  of  the  arc  O'P  onto  the  projection  plane. 
It  la  evident  froa  Figure  6  that 

(9)  O’P' .  2  R  tan  § 

To  reduce  equation  (9)  to  rectangular  coordinates.  on  the 
projection  plane,  the  coordinate  system  Illustrated  In  Figure  7 
Is  developed.  In  this  coordinate  system,  O'Y*  represents  the 
projection  of  the  meridian  through  O'  and  O'X*  represents  the 
projection . of  the  great  circle  through  O'  vhlch  Is  perpendicular 
to  the  meridian  through  O' . 

From  Figure  7,  the  following  Is  derived: 

X'  ■  F'P*  ■  O'P'  cos  B  ■  2  R  tan  ^  cos  B 

Y' •  I'P'»  O'P'  sin  B  ■  2  R  tan  §  sin  B 

(11) 

Utilising  the  spherical  triangle  formed  In  Figure  5, 
equations  (4),  (5),  and  (6)  can  again  he  derived.  By  substituting 
equations  (4),  (5),  and  (6)  Into  equations  (10'  and  (ll),  the 
mapping  equations  for  the  Stereographic  Projection  are  obtained: 

(12)  X'" 

c  c 


2  R  »in  q  sin  B 
!♦  cos  q 


a  2  R  sin  q  cos  B 
1  ♦  cos  q 
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(13) 


2  P  f cos  «  **n  f  -  sin  f  cos  f  cos  X) 

_ c _  c _ 

1  ♦  sin  |  'sin  #  ♦  cos  f  cos  f  cos  X 

c  c 

2.3  The  Orthographic  Projection 
The  Orthographic  Projection  is  developed  by  placing  a 
projection  center  at  infinity  and  by  projecting  the  datum  sphere 
onto  a  projection  plane  perpendicular  to  the  projection  rays.  All 
projection  rays  are  parallel  to  each  other  and  perpendicular  to  the 
projection  plane.  The  mapping  equations  for  the  Orthographic 
Projection  can  be  developed  in  analogous  manner  to  the  development 
of  the  mapping  equations  in  Section  2.1  and  Section  2.2.  Figure  8 
represents  a  sphere  with  a  radius  R  and  a  center  at  0.  The 
projection  center  is  located  at  infinity  in  the  direction  indicated 
by  the  figure.  The  projection  plane  is  perpendicular  to  the 
projecting  rays  and  tangent  to  the  sphere  at  O',  with  O'  also  being 
the  image  center  of  the  projection.  P  is  an  arbitrary  point  located 
on  the  sphere  at  a  latitude  f  and  at  a  longitude  difference  of  X 
from  O'.  P*  is  the  projection  of  P  on  the  projection  plane. 

Figure  8  illustrates  the  following: 

00*  -  OP  ■  R 
Angle  O'HP*  X 
Angle  RO'P*  90  -  B 
Arc  HO'  ■  90  -  # 

c 

Arc  RP  »  90  -  # 

Figure  9  i®  drawn  in  the  plane  of  triangle  OO'P  and 
illustrates  the  projection  of  the  arc  O'P  onto  the  projection  plane. 
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It  is  evident  from  Figure  9  that 

(lU)  O'P'-R  Bin  q 

To  reduce  equation  (l4)  to  rectangular  coordinates  on  the 
projection  plane,  the  coordinate  system  illustrated  in  Figure  10 
is  developed.  In  this  coordinate  system,  O'Y'  represents  the 
projection  of  the  meridian  through  O'  and  O'X*  represents  the 
projection  of  the  great  circle  through  O'  vhlch  is  perpendicular 
to  the  meridian  through  O' . 

From  Figure  10,  the  following  is  derived: 

(15)  X'aFP'a  O'P*  cos  B  ■  R  sin  q  cos  B 

(16)  Y'  -  EP*  -  O'P'  sin  B  -  R  sin'q  sin  B 

Utilizing  the  spherical  triangle  formed  in  Figure  8, 

equations  (4)  and  (5)  can  again  he  derived.  By  substituting 
equations  (U)  and  (5)  into  equations  (15)  and  (l6),  the  mapping 
equations  for  the  Orthographic  Projection  are  obtained: 

(17) 

(18) 


X'  »R  cos  f  sin  X 

Y'BR  (cos  |  sin  |  -  sin  #  cos  |  cos  X) 
c  c 
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3.  A  DESCRIPTIVE  GEOMETRY  APPROACH  TO  PERSPECTIVE  MAP  PROJECTIONS 

The  three  perspective  asp  projections  just  discussed,  the 
Qnoaonic,  the  Stereographic,  and  the  Orthographic,  were  developed 
by  placing  the  projection  center  at  a  particular  location  In 
reference  to  the  datum  sphere.  In  this  chapter  the  projection  center 
will  be  placed  at  an  arbitrary  distance  above  the  surface  of  the 
datua  sphere  and  perspective  aapping  equations  vill  be  developed  for 
this  spherical  datua  surface.  Since  the  projection  center  will  be 
located  at  an  arbitrary  distance  above  the  datum  surfo.ee.  It  Is  new 
possible  to  better  visualize  the  physical  significance  of  the 
perspective  projection  by  comparing  it  to  an  aerial  photograph.  As 
was  stated  earlier,  a  photograph  Is  actually  a  type  of  perspective 
projection.  Thus,  the  perspective  napping  equations  derived  In 
this  chapter  would  enable  the  plotting  of  grid  lines  on  an  Ideal 
vertical  photograph  taken  at  a  particular  distance  above  a  datum 
sphere.  Naturally,  this  can  not  be  accomplished  with  an  actual 
aerial  photograph  of  the  earth  due  to  the  fact  that  the  earth  Is 
not  a  sphere  nor  a  staple  geometric  figure  of  any  type.  Also,  due 
to  nuaerous  distortions  Inherent  to  photography,  aerial  photographs 
are  not  Ideal  vertical  projections.  However,  the  comparison  between 
perspective  projections  and  vertical  aerial  photographs  Is  still 
useful  in  visualising  the  derivations  Included  in  this  and  the 
follearlng  chapter. 
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Figure  11  illustrates  the  descriptive  geometry  approach  to 

perspective  map  projections  presented  in  [12].  In  this  figure  the 

lover  circle  represents  a  transverse  or  front  view  of  the  datum 

sphere,  with  the  plane  of  the  X  axis  defining  the  plane  of  the 

equator  and  the  Y  axis  defining  the  polar  axis.  The  upper  circle  in 

Figure  11  represents  the  polar  or  top  view  of  the  datum  Bphere, 

vlth  N  denoting  the  position  of  the  pole.  The  small  circle  on  the 

right  aide  of  the  figure  represents  the  perspective  view  of  the 

datum  sphere  that  would  be  obtained  by  placing  a  projection  center 

at  C  .  In  this  chapter  mapping  equations  will  be  derived  which 
F 

allow  the  transformation  of  coordinates  of  an  arbitrary  point 
from  a  three  dimensional  latitude -longitude  system,  as  illustrated 
in  the  top  and  front  views  of  Figure  11,  to  a  plane  X '  -Y*  system, 
as  illustrated  on  the  projection  plane  in  the  perspective  view  of 
Figure  11. 

In  Figure  11,  C  represents  a  projection  center  located  at 
F 

a  latitude  {  .  at  a  longitude  \  ,  and  at  a  distance  h  above  a 
c*  c 

datum  sphere  of  radius  R.  A  cene  ef  projection  rays  emanating  from  C_ 

is  tangent  to  the  sphere  at  A  and  B  and  creates  right  angles  O^AC^ 

and  C^BCy.  The  complete  intersection  of  the  projection  cone  and  the 

Bphere  is  a  circle  whose  trace  AB  is  shown  in  the  front  view.  This 

circular  intersection  as  seen  from  the  projection  center  C  is 

r 

illustrated  in  the  perspective  view.  This  circular  area  can  be 
visualized  as  that  portion  of  the  datum  sphere  which  could  be  seen 
if  an  observer  were  located  at  C  .  In  this  derivation  the  AB 
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plane  will  be  used  as  the  projection  plane  and  all  points  on  the 

visible  portion  of  the  sphere  will  be  projected  onto  this  plane. 

In  the  front  view  of  Figure  11,  iet  P  be  an  arbitrary 

point  located  at  a  latitude  f  and  at  a  longitude  X  on  the  portion 

of  the  sphere  visible  from  Cp.  By  definition,  X  equals  the 

difference  in  longitude  between  the  arbitrary  point  P  and  the 

projection  center  C.  P^  represents  the  position  of  the  point  P  in 

the  top  view.  P  is  projected  along  the  projection  ray  C  P  until 
*  F  F 

the  ray  intersects  the  projection  plane  A  at  P_‘.  P  '  is  then 

'T  F 

the  projected  image  of  Pp  on  the  projection  plane.  P  '  represents 

the  position  of  this  point  in  the  top  view.  It  is  seen  that  the 

distance  DPp'  is  then  the  Y*  coordinate  of  the  point  P'  on  the 

projection  plane.  In  the  top  view,  the  distance  P  '  *P  '  is  the  X' 

T  T 

coordinate  of  the  point  P'  on  the  projection  plane  and  this  distance, 
when  transferred  to  the  perspective  view,  locates  the  point  P'. 

P',  with  an  X*  value  of  P  ’ 'P  ’  and  a  Y'  value  of  DP  ' ,  is  thus 

the  projection  of  point  P  on  the  projection  plane  in  the  perspective 

view. 


The  point  P  has  now  been  transferred  to  the  projection  plane, 
and  it  only  remains  necessary  to  resolve  the  ordinate  and  abscissa 
values  of  the  point  into  terms  involving  the  original  given  para¬ 
meters.  This  is  accomplished  in  the  following  paragraphs. 


The  angle  C  AO  is  designated  0  and  it  is  evident  that 
F  f 


(19) 


cos  6 


AO 
_ F 

°FCF 


R 

R  ♦  h 
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I 


and  also  that 

(20)  0_D ■  R  cose 

r 

The  equation  of  the  line  AB  may  now  he  written  as 

(21)  X  cos  #  +Y  sin  I  -  R  cos  Q  ■  0 

c  c 

and  the  equation  of  the  line  may  he  written  as 

(22)  X  sin  (  -  Y  cos  t  ■  0 

c  c 

Next,  the  two  perpendicular  distances  LP  and  P  M  from 

F  F 

point  P_  to  lines  AB  and  0  C  respectively  must  he  determined.  This 
F  F  F 

is  done  hy  using  well  known  geometric  formulas  for  finding  the 
perpendicular  distance  from  a  point  to  a  line.  The  resulting 
equations  are 

(23)  LP  z  R  (c<>8  4  cob  t  cos  X  +  sin  |  sin  f  -  cos  e) 

F  c  c 

(24)  MP  ■  -R  (sin  #  cos  |  cos  X  >ces  f  sin  |) 

F  c  c 

Considering  the  X^Y*  coordinate  system  of  the  projection 

plane  in  the  perspective  view,  it  was  shown  previously  that  the 

ordinate  of  P*  in  this  coordinate  system  was  equal  to  DP '  ' , 

F 

Utilizing  similar  triangles  in  the  front  view  the  fallowing  is 
attained: 

MP  DC  MP 
F  F  F 

DP  '  »DC  - fp — 

F  F  "'f  ^f  ’  “V 

and  also  from  the  front  view  the  following  is  derived: 

DC  «AD  tan  0  ■  R  sin  0  tan  0 
F 

A  combination  of  the  above  equations  yields  the  Y*  cordinate 
of  P'  in  the  X’-Y'  coordinate  system  of  the  projection  plane. 
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1 


(R  sin  8  tan  0)  R  (cos  #  sin  |  -  sin  #  cos  |  cocT  >.) 
R  sin  0  tan 0  -  R(cos  f  cos  f  cos  A+sln  f  sin  #-  cos  g) 


(25) 


2 

R  sin  6  (cos  f  sin  #  -  sin  f  cos  f  cos  \) 


i-  c  c 

Y '  1  -  cose  (cos  •  cos  f  cos  X+  sin  ?  sin  #) 

c  c 


A  X'  coordinate  corresponding  tc  the  Y’  coordinate  calculated 
by  equation  (25)  may  now  be  derived.  Considering  Figure  11  again, 
the  following  proportions  are  obtained: 


P  up  »  p  *C  p  »C  DP  • 
T  T_T  T  _  F  F  _  F 

p--'g  'rr  rr  r 

T  T  t  T  F  F  F 


It  is  also  evident  from  Figure  11  that 


"R  cos  «  sin  \ 

T  T 


DPy '  •  Y ' 


The  abscissa  X'  of  the  point  P'  on  the  projection  plane 
was  previously  shown  to  be  equal  to  P^’P^'  ’ind  thu8 


p  i.-p  DP  ' 
<p  T  *  F 


X«-  p  • ip  i  ■ 

T  T  HP_ 


■  R  cos  ♦  sin  \  »  Y* 

X'  -R  (sin  |  cos  ff  cos  sin  #  sin  #) 


.  ■ _  R  sin  8  cos  t  sin  X 

x  1  -  cos  6  (cos  |  cos  f  cos  X ♦  sin  |  sin  |) 


Summarizing,  equations  (25)  and  (26)  are  perspective 


projection  sapping  equations  which  enable  an  arbitrary  point  to  be 
transferred  from  a  three  dimensional  latitude-longitude  coordinate 
system  of  a  datum  sphere  to  a  plane  X'-Y'  coordinate  system  of  a 
projection  plane. 
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4.  A  PHOT OGRAEtffiTR IC  APPROACH  TO  PERSPECTIVE  MAP  PROJBCTIOWS 


Equations  (25)  and  (26),  developed  by  a  descriptive  geaoetry 
method  in  Chapter  3>  are  mapping  equations  suitable  for  a  spherical 
datum  surface.  It  is  not  apparent  to  the  author  hour  these  equations 
can  be  modified  to  make  them  applicable  to  an  ellipsoidal  datum 
surface.  It  is  also  not  immediately  apparent  how  these  equations 
can  be  used  for  all  possible  locations  of  the  projection  center.  In 
this  chapter  mapping  equations  will  be  derived  that  can  be  utilised 
regardless  of  the  location  of  the  projection  center  and  that  can  be 
used  both  on  spherical  and  ellipsoidal  datum  surfaces. 

Since  a  perspective  projection  is  directly  comparable  to  a 
photograph,  the  author  chose  to  derive  the  desired  mapping  equations 
through  a  photogrammetric  approach.  Section  4.1  deals  with  the 
derivation  of  equations  based  on  a  spherical  datum  surface,  while 
Section  4.2  deals  with  similar  equations  based  on  an  ellipsoidal 
datum  surface.  Since  a  sphere  is  only  a  special  case  of  an 
ellipsoid,  napping  equations  derived  for  the  sphere  are  only  special 
cases  of  those  derived  for  the  ellipsoid.  The  author  considers 
the  spherical  case  worthy  of  separate  discussion  based  on  the  belief 
that  the  sphere  Is  a  logical  intermediate  step  before  progressing 
to  the  more  complicated  ellipsoid.  The  mapping  equations  derived 
for  the  sphere  may  also  be  compared  with  previously  derived  equations 
as  a  check  on  the  validity  of  the  general  method  being  utilized  in 
this  chapter. 
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4.1  A  Spherical  Datum  Surface 
4.11  The  Mapping  Equations 

Figure  12  represents  a  spherical  datum  surface  vith  a  radius 
R  and  a  conventional  three  dimensional  X-Y-Z  coordinate  system. 
Point  C  depicts  the  projection  center  and  Is  located  at  a  latitude 
#c,  at  a  longitude  Xc,  and  at  a  height  above  the  datum  sphere  of 
h.  Point  P  represents  an  arbitrary  point  on  the  datum  surface 
and  Is  located  at  a  latitude  #  and  at  a  longitude  X.  In  the 
following  derivation  X  will  equal  the  difference  in  longitude  between 
the  point  P  and  the  projection  center  C. 

The  space  rectangular  coordinates  of  point  C  are  obtained 
from  Figure  12. 

(27)  Y  .  (R*h)  cos  f 

*  C 

(28)  xc.0 

(29)  2c,(R.h)  .in  #c 

Likewise,  the  coordinates  of  point  P  are  obtained  from  Figure  12. 

(30)  Y  »R  cos  |  cos  X 

P 

o 

(31)  X  ■  R  cos  f  sin  X 

P 

(32)  Z-Rsinf 

P 

As  an  aid  In  visualizing  the  remainder  of  the  derivation, 
point  C  may  be  compared  to  a  camera  lens  nnd  CO,  the  projection 
axis,  Is  comparable  to  a  camera  axis.  The  projection  plane  AB, 
or  mapping  plane,  may  be  visualized  as  a  photographic  plate  which 
is  perpendicular  to  the  camera  axis  and  located  at  a  distance  f 
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Figure  13 
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behind  the  lens.  It  is  seen  that  f  is  comparable  to  the  fecal 
length  of  a  camera  and,  thus,  the  scale  of  the  projection  is 
dependent  on  f .  Actually,  the  ratic  of  f  to  h  determine*  the 
•cale  of  the  projection;  and  it  is  apparent  that  a  full  scale 
projection  will  result  vhen  f  is  equal  to  h.  It  will  be  necessary 
later  in  the  derivation  to  have  firm  definitions  of  both  h  and  f 
as  regards  to  sign  convention.  The  signs  of  h  and  f  are  defined 
as  positive  vhen  they  are  above,  or  to  the  outside  of,  that  portion 
of  the  sphere  to  be  projected  and  negative  if  to  the  inside  of  the 
sphere. 

Figure  12  illustrates  the  following: 

CS-h 

CC'»  f 

06  -  OP-  R 

The  next  step  in  the  derivation  Is  to  rotate  the  original 
Z  axis  in  the  Z-Y  plane  until  the  Z  axis  coincides  vlth  the 
projection  axis  0C,  as  is  illustrated  in  Figure  13.  This  rotation 
in  the  Z-Y  plane  creates  a  new  X'-Y'-Z'  coordinate  system.  It  Is 
evident  that  this  type  of  rotation  will  not  alter  the  X  coordinate 
of  points  from  their  values  as  calculated  in  the  original  coordinate 
system. 

Figure  13  illustrates  the  follewing: 

PL«  Y 
P 

PU.  Z 

P 
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PK-Y  • 

P 

PM-Z  • 

P 

NO-Y  sin  « 
P  c 


UT-  Z  cos  f 
P  c 

UN-  Y  cos  # 
P  c 


PT-Z  sin  | 
P  c 


Considering  Figure  13 ,  the  space  rectangular  coordinates  of 


point  C 

in  the  new  X'-Y'-Z'  coordinate  system  are  obtained. 

(33) 

X  •  -  0 
c 

(34) 

Y  '  -0 
c 

(35) 

Z  '  ■  R+  h 
c 

Likewise, 

the  new  coordinates  of  point  P  are  obtained. 

(36) 

Y  '  ■  Y  sin  #  -  Z  cos  f 

P  p  c  p  c 

(37) 

X  X 

P  P 

(38) 

Z  '  ■  Y  cos  |  +Z  sin  f 

P  p  c  p  c 

Substituting  equations  (30),  (31),  and  (32)  into  equations  (36),  (37), 
and  (38),  the  X'-Y'-Z'  coordinates  of  point  F  in  terms  of  the  original 
parameters  are  obtained, 

(39)  Y  '  ■  R  sin  #c  cos  |  cos  \  -  R  cos  #c  sin  | 


27 


• 

(40) 

X  •  - 
p 

R  cos 

4  sin  \ 

f 

(M) 

z  •• 

p 

R  cos 

f  cos  4  c os  UR  sin  4  sin  f 
c  c 

The  equation  of  the  projection  plane  AB  which  is  perpendicular 
to  the  projection  axis  and  hence  now  perpendicular  to  the  Z’  axis 
is  clearly 

(42)  Z»  .  R  ♦  h  ♦  f 

The  equation  of  the  projection  ray  from  P  that  passes 
through  C  and  intersects  the  projection  plane  at  P'  can  be  written 

A8 


Equation  (43)  is  the  equation  of  the  line  in  space  along 
which  P  is  projected  onto  the  projection  plane.  The  intersection 
of  this  line  and  the  projection  plane,  as  given  by  equation  (42), 
is  the  location  of  P',  the  projection  of  P  o*  the  projection  plane. 
Therefore,  the  simultaneous  solution  of  equations  (42)  and  (43) 
will  yield  the  desired  X'-Y'  coordinates  of  the  projected  point. 

The  equation  for  X*  is  derived  in  the  following  manner* 
From  equation  (43) 

X*  -  X  •  Z»  -  Z  * 
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m 


x»- 


(Z*  -  Z  •  X  •  -X  •) 

- iJ  i  \ — E— +x  • 

(V  -  z  »)  p 


(45) 


Substituting  equation  (33)  into  equation  (44) 
X' 


V  (V  -z,> 


-C  p 

Substituting  equations  (35),  (4o),  (4l)  and  (42) 
into  equation  (45) 

-  f  R  cos  t  sin  X 

x'“  R  (1  -  cos  *  cos  f  cos  X  -  sin  #  sin  #)  ♦  h 
c  c 

Remembering  that  a  photogrammetric  approach  was  utilised 
in  this  derivation,  it  will  be  realized  that  the  X'  coordinate  has 
suffered  a  reversal  of  sign  during  the  process,  presenting  a  reversed 
image  of  the  original  datum  surface.  In  order  to  compensate  for  this 
situation  the  sign  of  the  above  equation  must  be  reversed,  and  there¬ 
fore 


(46) 


X'  . 


f  R  ccs  |  sin  X 


_  *>  *  ■  V'-W  V  *  - _ 

R  (1  -  cos  f  cos  f  cos  \  -  sin  #  sin  |)  +  h 
c  c 


The  equation  for  Y‘  is  derived  in  a  manner  analogous  to  the 
derivation  of  the  equation  for  X’  as  follows: 

From  equation  (fcj) 


(47) 


Y»  -  Y  »  Z»  -  Z  « 

- _ Z-l£ — 

Y  •  -  Y  *  Z  '  -  Z  ' 
c  p  c  p 


(z-  .  Z  •)  (Y  -  -  Y  •) 

11  ■  — (V  -X’) p  *  h 
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Substituting  equation  (3*0  into  equation  (47) 


(48) 


YS.  V  <V  -Z’> 

*e.  -  V 


Substituting  equations  (35),  (39),  (4l),  and  (42)  into 

equation  <48) 


f  R  (cos  fc  sin  |  -  sin  #c  cos  f  cos  X) 

^9)  R  (l  -  cos  *  cos  #  cos  X  -  sin  $  sin  #)  ♦  h 

c  c 

In  summary,  equations  '46)  and  (49)  are  perspective 
projection  mapping  equations  derived  for  a  spherical  datum  surface. 
These  equations  allow  the  calculation  of  plane  X'-Y'  coordinates  of 
a  point  P'  on  a  projection  plane  located  at  a  distance  f  from 
the  projection  center.  Initially,  point  P  had  been  located  on  the 
datum  sphere  by  spherical  #-X  coordinates  and  C  by  both  spherical 
#  -X  coordinates  and  a  distance  h  above  the  datum  surface. 

C  w 


4.12  A  Comparison  Of  Mapping  Equations 

Two  sets  of  mapping  equations  have  now  been  developed  for 
the  case  of  a  projection  center  being  placed  at  an  arbitrary 
distance  above  a  spherical  datum  surface.  If  both  sets  of 
equations  are  valid,  they  must  prove  to  be  equal.  A  brief 
comparison  of  equations  (25)  and  (26),  derived  by  the  descriptive 
geometry  approach  of  Chapter  3,  with  equations  (46)  and  (49),  Just 
derived,  illustrates  a  great  similarity  in  the  equations.  The 
differences  in  the  equations  are  due  to  the  fact  that  while 
equations  (46)  and  (49)  contain  the  parameters  h  and  f,  equations 
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(25)  and  (26)  contain  the  comparable  parameter  q.  In  effect,  the 
parameter  9  fixes  a  f  for  a  particular  h.  As  explained  earlier, 
however,  since  f  acts  only  as  a  scale  factor,  its  value  is  not  of 
real  significance. 

Considering  Figure  11  and  recalling  that  AB  was  the  pro¬ 
jection  plane,  it  is  evident  that  DC,  being  the  distance  between 
the  projection  center  and  the  projection  plane,  is  the  parameter 
f  introduced  in  Section  4.11.  From  Figure  11,  the  following 
expression  for  the  parameter  f  may  be  derived: 

f "  "^F^F  "  ®F^ 

■  R  +  h  -  R  cos  9 

2 

(50)  -  R  +  h  -  R  +  h 

Utilizing  the  definition  of  the  parameter  9  as  given  by 
equation  (19)  and  the  expression  for  f  as  given  by  equation  (50), 
the  perspective  projection  mapping  equations  derived  by  the 
descriptive  geometry  method  of  Chapter  3  will  be  proved  equal  to 
the  mapping  equations  derived  by  the  photogrammetric  method  of 
Section  4.11.  Equation  (26),  derived  in  Chapter  3,  nay  he  trans¬ 
formed  into  equation  (46),  derived  in  Section  4.11,  as  follows: 

Rewriting  equation  (26) 

2 

_ R  sin  6  co»  ♦  sin  A 

x' *  1  -  cos  6  (cos  |  cos  f  cos  A ♦  sin  f  sin  |) 

c  c 

2  2 

Substituting* (l  -  cos  ©)  for  sin  q  and  equation  (19) 
into  the  above  equation 
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[R ♦  h-  R'  "r —  ]R  cos  y  sin  X 

V  y  *— m~w — • _  - - -  —  - - -  _  - ^ 

R  (l  -  cos  |  cos  |  cos  X.  -  sin  y  sin  $)  *  h 
c  c 

Substituting  equation  (50)  into  the  above  equation 


(51) 


X’. 


f  R  cos  f  sin  X 


*  *»  v  vw  x  w  n 

R  (l  -  cos  y  cos  #  cos  X  -  sin  ft  sin  y)  ♦  h 


In  a  like  manner,  equation  (25),  derived  in  Chapter  3>  roay  be 

transformed  into  equation  (49),  derived  in  Section  4.11,  as  follows: 

Rewriting  equation  (25) 

2  . 

R  sin  6  (cos  y  sin  y  -  sin  f  cos  y  cos  X) 
c  c 

Y'“  l  -  cose  (cos  y  cos  y  cos  X  +  sin  y  sin  y) 

e  o 


2  ? 

Substituting  (l  -  cos  0)  for  sin  ©and  equation 


(19)  into  the  above  equation 


d 

(R  +  h  -  _JL —  )  R  (cos  y  sin  y  -  sin  y  cos  y  cos  X) 

Y t  .  .I,  R  t  ^  1 .  ..  S  .9 ... 

*  R  (1  -  cos  y  cos  y  cos  X  -  sin  y  sin  y)  +  h 
c  f 


Substituting  equation  (50)  into  the  above  equation 


f  R  (cos  yc  sin  y  -  sin  yc  cos  y  cos  X) 

(52)  “  R  (l  -  cos  y  cos  y  cos  X  -sin  y  sin  y)  +  h 

c  c 

Equations  (51)  Rnd  (52)  are  seen  to  be  identical  to 
equations  (46)  and  (49).  Thus,  the  perspective  nrojection  mapping 
equations  derived  by  the  descriptive  geometry  method  are, 
in  effect,  equal  to  the  equations  derived  by  the  photogrammetric 
method.  This  equality  is  a  partial  proof  of  the  validity  of  the 
derivations  and  the  derived  equations  of  Section  4.11. 


32 


4.13  The  Gncmonic  Projection 

Three  special  cases  of  the  perspective  projection  were 
discussed  in  Chapter  2,  and  three  sets  of  mapping  equations  were 
developed  for  these  special  cases.  Since  the  perspective  projection 
mapping  equations,  equations  (46)  and  (49),  are  applicable  to  all 
perspective  projections  of  a  spherical  datum  surface,  these  equations 
must  be  applicable  also  to  the  above  mentioned  special  cases.  This 
fact  will  be  proved  in  this  and  the  two  succeeding  sections. 

As  was  illustrated  in  Section  2.1,  the.Gnomonic  Projection 
is  a  special  case  of  the  perspective  projection  in  which  the 
projection  center  is  located  at  the  center  of  the  sphere.  By 
reviewing  Figure  2  and  by  recalling  the  definitions  given  earlier 
for  h  and  f,  the  values  for  h  and  f  in  the  Gnomonic  Projection  are 
seen  to  be  -R  and  -R  respectively.  Substituting  these  values  into 
the  perspective  projection  mapping  equations,  equations  (46)  and 
(49),  yields  the  following: 

-R  R  cos  $  sin  \ 

X'"  R  (l  -  cos  i  cos  $  cos  X  -  sin  $  sin  #)  -  R 
c  c 

•  K  cos  |  sin  A 

cos  *  cos  i  cos  \  +  sin  ft  sin  f 
c  c 

“  R  R  (cos  |  sin  f  -  sin  f  cos  $  cos  A) 

_ c _ c _ _ _ 

"  R  (l  “  cos  f  cos  f  cos  A  •*  sin  f  sin  #)  -  R 
c  c 

R  (cos  ft  sin  f  -  sin  #  cos  |  cos  A) 

m  C  C 

cos  *  cos  f  cos  \  ♦  sin  f  sin  # 
c  c 


(53) 


(54) 
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Equations  (53)  and  (54)  are  identical  to  equations  (7)  and 
(8)  developed  for  the  Gnomenic  Projection  In  Section  2.1.  Thus,  the 
napping  equations  for  the  Gnooonic  Projection  have  been  proved  to 
be  special  cases  of  the  perspective  projection  mapping  equations. 

4.14  The  Stereographic  Projection 

The  Sterographic  Projection  was  the  second  special  case  of 
the  perspective  projection  discussed  in  Chapter  2.  Considering 
Figure  5 r  the  values  for  h  and  f  in  tu«'  Stereographic  Projection 
are  seen  to  be  -2R  and  -2R  respectively.  Substituting  these  values 
into  the  perspective  projection  mapping  equations,  equations  (46) 
and  (49),  yields  the  following: 

-2R  R  cos  |  sin  X _ 

x "  ft  (l  -  cos  cos  #  cos  k  -  sin  |  sin  |)  -  2R 
c  c 

■  2R  cos  ft  sin  1 

1  +  cos  |  cos  f  cos  X  ♦  sin  #  sin  * 
c  c 

-2R  R  (cos  fc  sin  |  -  sin  #c  cos  «  cos  X) 

R  (1  -  cos  f  cos  f  cos  X  ♦  sin  #  sin  #)  -  2R 
c  c 

2R  (cos  |  sin  •  -  sin  #  cos  |  cos  X) 
s  ■ 

1  ♦  cos  lc  cos  |  cos  X  ♦  sin  9C  sin  f 

Equations  (55)  and  (56)  are  identical  to  equations  (12)  and 
(13)  developed  for  the  Stereographic  Projection  in  Section  2.2. 

Thus,  the  equations  for  the  Stereographic  Projection  have  been  proved 
to  be  special  cases  of  the  perspective  projection  mapping  equations. 


(55) 


(56) 


U.15  The  Orthographic  Projection 

The  third  special  case  of  the  perspective  projection 
discussed  in  Chapter  2  was  the  Orthographic  Projection.  In  order 
to  convert  the  perspective  projection  mapping  equations  for  use  in  an 
Orthographic  Projection,  a  relationship  must  he  established  between 
f  and  h.  It  was  discussed  earlier  that  f  has  the  effect  of  a 
scale  factor,  and  from  elementary  projection  principles  it  was 
seen  that  a  full  scale  projection  results  if  f  equals  h.  Therefore, 
in  developing  an  Orthographic  Projection,  f  may  be  set  equal  to 
h.  Making  this  substitution  in  the  perspective  projection  mapping 
equations,  eqaations  (U6)  and  (^9),  yields  the  following: 


The  numerators  and  denominators  of  the  above  equations  are  now 

t 

divided  by  h,  producing: 


R  coo  4  sin  X 

R  (1  -  ccs  4^  cos  4  co3  X  -  sin  4„  sin  #' 


R  (cos  4^  sin  4  -  sin  4fi  cos  4  cos  X) 


R  vl  -  cos  4  cos  4  cos 


sin  4 sin  4J 
c 


The  Orthographic  Projection  is  based  on  the  projection  center  being 
located  at  infinity.  Substituting  infinity  for  h  in  the  above 
equations  yields: 
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(57)  X*  ■  R  cos  f  sin  X 

(58)  Y’  ■  R  (cob  f  sin  f  -  sin  |  cot  I  cos  X) 

c  5 

Equations  (57)  nnd  (58')  are  Identical  to  equations  (IT)  and 
(l8)  developed  for  the  Orthographic  Projection  in  Section  2.3. 

Thus,  the  equations  for  the  Orthographic  Projection  have  also  been 
proved  to  be  special  cases  of  the  perspective  projection  mapping 
equations. 


4.2  An  Ellipsoidal  Datum  Surface 
4.21  The  General  Mapping  Equations 

Section  4.1  dealt  vith  perspective  projection  mapping 
equations  for  a  spherical  datum  surface;  however,  the  sphere  is 
actually  a  rather  poor  approximation  for  the  shape  of  the  earth.  A 
better  approximation  for  the  shape  of  the  earth  is  an  ellipsoid  of 
revolution,  and  this  section  deals  with  the  derivation  of  General 
Perspective  Projection  Mapping  Equations  Tor  ellipsoidal  datum 
surfaces.  In  the  remainder  of  this  paper,  an  ellipsoid  of  revolution 

will  be  implied  by  the  term  ellipsoid.  The  parameters  used  to  define 

2 

the  ellipsoid  will  be  a,  the  major  semi -diameter,  and  e  ,  the  square 
of  the  eccentricity.  It  will  be  seen  that  the  derivation  method 
to  be  used  in  this  section  is  analogous  to  the  method  employed  in 
Section  4.1L 

Figure  l4  illustrates  the  datum  surface  which  is  an  ellipsoid 
of  revolution  with  the  Z  axis  being  the  axis  of  rotation.  Point  C 
depicts  the  projection  center  and  is  located  at  a  geodetic  latitude 


36 


p* 

cy 


z 
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cp  ,  at  a  longitude  X  ,  and  at  a  height  above  the  datum  surface  of 
h.  Point  P  represents  an  arbitrary  point  on /the  datum  surface  and 
is  loc  ted  at  a  geodetic  latitude  cp  and  at  a  longitude  X.  As  before, 
X  will  equal  the  difference  in  longitude  between  the  point  P  and 
the  projection  center.  Nc  a.id  Np  are  the  radii  ol  curvature  in  the 
prime  vertical  at  C  and  P  respectively. 

Figure  14  illustrates  the  following: 

CS»  h 
CC'.  f 
KP  -  Np 

O'S-  N 

c 

The  space  rectangular  coordinates  of  point  C  are  obtained 
from  Figure  14. 

(59)  Xc-o 

(60)  Y  ■  (N  ♦  h)  cos  cp 

C  C  c 

(61)  Z  •  [N  (l  -  e  )  ♦  h)  sin  cp 

Likewise,  the  coordinates  of  point  P  are  obtained  from  Figure  lb. 

(62)  X  •  N  ^os  cp  sin  X 

p  P  y 

(63)  Y  ■  N  cos  cp  cos  X 

P  P 

(64)  Z  •  N  (1  -  e2)  sin  m 

P  P 

Similar  to  the  derivation  utilizing  a  spherical  datum 
in  Section  4.11,  the  Z  axis  is  now  rotated,  and  in  this  case 
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transposed,  in  the  Z-Y  plane  until  the  Z  axis  coincides  with  the 
projection  axis  CO'.  The  new  X'-Y'-Z'  coordinate  system  is 
illustrated  in  Figure  15.  The  values  of  X  coordinates  of  points 
are  not  changed  by  the  rotation  due  to  the  rotation  being  accomplished 
in  the  Z-Y  plane. 

Figure  15  illustrates  the  following: 


PH  *Y 

P 

PU-Z 

P 

PL-Y  • 

P 

PM-Z  ' 

P 

JO^Y  sin 

fP 

P 

c 

TK*=  00'  cos 

*c 

UR*  Z  cos 

P 

c 

JT  «■  Y  cos 

<p 

P 

c 

UK  -  00'  sin 

RPkZ  sin 

9 

P 

c 

The  distance  00’  which  the  Z'  axis  has  been  transposed  along 

the  Z  axis  is  computed  ns  follows: 

00' c  (N.  ♦  h)  sin  u  -  Z 
c  Tc  c 

■(Nc  ♦  h)  sin  9c  -  [Nfi  (l  -  e^)  +  h]  sin  cp^ 

(65)  ’  Hc  ein  % 
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The  space  rectangular  coordinates  of  C  in  the  new  X'  -Y*  -Z' 
coordinate  systen  are  obtained  from  Figure  15. 


(65®) 

V*° 

(66) 

m 

■ 

o 

(67) 

Z  •  •  H  ♦  h 
c  c 

Likewise,  the  new  coordinates  of 

point  P  are  obtained. 

(68) 

X  •  -  X 

P  P 

(69) 

Y  »  -  Y  sin  cp  -  (Z  ♦  00») 
P  p  c  p 

cos  cp 
c 

(70) 

Z  •  •  Y  cos  cp  ♦  (Z  ♦  00' ) 
P  P  c  p 

sin  cp 
c 

Substituting  equations  (62),  (63),  (6U), 

end  (65)  into  equations 

(68), 

(69),  and  (70),  the  coordinates  of 

point  P  are  obtained 

in  terms  of  the  original  parameters. 

(71) 

X  *  ■  H  cos  cp  sin  X 

P  P 

Y  »  ■  H  sin  cp  cos  cp  cos  X 
P  p  c 

,  2. 

-  H  (1  -  e  )  cos  cp  sin 
P  c 

o 

(72)  -He  sin  cp  cos  cp 

c  c  c 


Z  •  »N  job  cp  cos  cp  cos  X+F  (l  -  e  )  sin  cp  sin  cp 
p  p  Tc  T  p  c 

(73)  **c  e2  sin2  9C 

The  equation  of  the  projection  plane  AB  which  is  perpendicular 
to  the  projection  axis  and  hence  now  perpendicular  to  the  Z*  axis 


is  now  foraulated  as 


The  equation  of  the  projection  ray  from  P  that  passes 
through  C  and  intersects  the  projection  plane  f,t  P’  is  now 
written  as 


(75) 


X«  -  X  •  Y*  -  Y  • 

•  -  fc'  *Y  '  -  Y  *  “Z  '  -  I  • 


"c  ”p  ~c  Yp  C 

Similar  to  the  derivation  in  Section  4.11,  the  simultaneous 
solution  of  equations  (75)  and  (74)  vil}  yield  the  desired  X'-Y' 
coordinates  of  the  projected  point  on  the  projection  plane. 

The  equation  for  X*  is  derived  in  the  following  manner; 

From  equation  (75) 


Z»  -  Z  • 

b 

P 


(76) 


(77) 


(Z*  -  Z  •)  (X  •  -  X  •) 

x'- - -  P  .l. — E_  ♦x  • 

(Z  '  -  z  •)  p 

c  p 

Substituting  equation  (65a)  into  equation  (76) 


X  «  (Z  •  -  Z') 

Yt.  p  c 

z  '  -  Z  ' 


■c  -p 

Substituting  equations  (67),  (71),  (73),  and  (74) 
into  equation  (77) 

-  f  N  cos  rp  sin  X 

X'»  — g - g  .  P...  .1- 

Nc  (l  -  e  sin  cpc)  -  Np  [cos  cos  9  coo  X 


—  z- 

+  (1  -  e  )  sin  ^  sin  q  ]  ♦  h 

As  in  Section  4.11,  the  sign  of  the  above  equation 
must  be  reversed  to  compensate  for  the  inversion  due  to  the 

photogrammetric  derivation,  and  therefore 

4l 


The  equation  for  Y'  is  derived  in  a  manner  "mlogous  to  the 


derivation  of  the  equation  for  X'  as  follows: 
Prom  equation  (75) 


Substituting  equation  (66)  into  equation  (79) 

(8o) 


Substituting  equations  (67),  (72),  (73),  *nd  (Jk)  into 
equation  (00) 


(01) 


-  sin  cpc  cos  9  cos  \  j  ] 

A 

♦  (l  -  e  )  sin  9C  sin  cp]  i  h 


In  summary ^equations  (7C>)  md  (8l)  ire  General  Perspective 
Projection  Mapping  Equations  derived  for  an  ellipsoidal  datura 
surface.  These  equations  allow  the  calculation  of  plane  X*-Y' 


coordinates  of  a  point  P»  on  a  projection  plane  located  at  a  distance 


f  from  the  projection  center.  Initially,  point  P  had  been  located 

on  the  datum  ellipsoid  by  geodetic  cp-X  coordinates  and  C  by  both 

geodetic  9  -X  coordinate  and  a  distance  h  above  the  datum  surface, 
c  c 

4,22  A  Comparison  of  Mapping  Equations  For 
The  Ellipsoid  And  For  The  Sphere 

A  sphere  is  actually  a  special  case  of  an  ellipsoid  in 
vhlch  the  ellipsoid  flattening  is  equal  to  zero.  Since  equations 
(78)  and  (8l)  were  derived  as  General  Perspective  Projection 
Mapping  Equations  for  an  ellipsoidal  datum  surface,  these  equations 
should  prove  to  be  suitable  for  a  spherical  datum  surface.  There¬ 
fore,  equations  (46)  and  (49),  derived  for  a  spherical  datum 
surface,  sheuld  prove  te  be  special  cases  a?  equations  (78)  and 
(8l).  This  is  Illustrated  as  fellows* 

Rewriting  equations  (78)  and  (8l) 

XVI 

*  (1  -  e*)  Sin  •  sin  pj+h 

C 

2  2 
f  toe  sin  9  ces  9  ♦  H  [(l  -  e  )  ces  9  sin  9 

c  C  C  T>.  e 

2  2 

Ec  (l  -  e  sin  9c)  -  [ces  9e  ces  9  cos  X 
-  sin  9  ces  9  cos  X  ]  ] 

“7 — \ ~ — 

♦  (1  -  e  )  sin  9c  sin  ^  ♦  h 

*3 


Considering  the  parameters  in  the  above  equations,  it  is  apparent 
the  following  is  true  when  the  datum,  surface  is  spherical: 

N  -  N  -R 
c  p 

e  •  0 

The  substitution  of  the  above  equalities  into  equations  (78)  and 
(81)  transforms  these  equations  into  equations  (46)  and  (49),  the 
sapping  equations  for  a  spherical  datum.  Thus,  equations  (78) 
and  (81)  are  valid  both  for  spherical  and  ellipsoidal  datum, 
surfaces.  Further,  since  it  was  shown  that  the  mapping  equations 
for  the  Gnomonic,  Stereographic,  and  Orthographic  Projections 
were  only  special  eases  of  the  mapping  equations  for  a  spherical 
datua,  then  they  too  are  only  special  cases  of  equaticnr  (^S)  and 
(8l).  Therefore,  equations  (78)  and  (8l)  are  truly  General 
Perspective  Projection  Mapping  Equations. 
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*A"  PROJECTION  DISTORTIONS 


Map  distortions  is  an  extensive  subject  and  no  attempt  will 
be  made  in  this  paper  to  treat  the  subject  in  general.  The  author 
has  derived  an  empirical  method  to  calculate  t.he  linear  distortions 
of  the  meridians  and  parallels  on  the  projection  plane  nd  the 
angular  distortions  of  the  intersections  of  the  meridians  and 
parallels  on  the  projection  plane.  The  calculation  of  these 
distortions  not  only  enables  a  vizualization  of  the  projected 
latitude-longitude  grid  but  also  may  be  used  as  indicators  ns  to 
the  type  of  distortions  that  will  be  found  in  projecting  *ny  line 
from  the  datum  surface  onto  the  projection  plane. 

5.1  Meridian  And  Parallel  Linear  Distortions 
The  author  defines  a  linear  distortion  to  exist  when  the 
linear  distance  between  two  points  measured  on  the  d-’tum  surfac 
does  not  equal  the  linear  distance  between  the  same  two  points 
projected  full  scale  onto  the  projection  plane.  The  amount  of 
distortion  will  be  indicated  by  a  ratio  of  the  distance  measured  on 
the  projection  plane  to  the  distance  measured  on  the  datum  surface. 
In  order  to  deterrine  the  linear  distortions  at  a  point,  the  two 
distances  to  be  compared  would  of  necessity  have  to  be  infinitely 
small.  However,  as  will  be  shown  by  examples  in  Chapter  6,  good 
approximations  can  be  made  by  utilizing  finite  distances. 
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The  distance  between  two  point*,  P.,  and  Pg,  located  by 
plane  X'-Y*  coordinates  on  a  projection  plane  is  given  by  the 
following  elementary  equation: 

(82)  d-  [(X2*  -  Xi')2*  (Yg*  -  Y1*)2]  * 

The  dlatance  between  two  points  on  the  sane  parallel  on  the 
projection  plane,  d^',  would  thus  be  obtained  by  using  equations 
(78)  and  (8l)  to  compute  the  X'-Y'  coordinates  of  the  two  points 
and  then  by  using  equation  (82)  to  calculate  the  distance  between 
the  points.  Likewise,  the  distance  between  two  points  on  the  same 

meridian  on  the  projection  plane,  d  ' ,  would  also  be  obtained  by 

m 

using  equations  (78),  (8l),  and  (82). 

The  distance  between  two  points  on  a  Meridian  on  the 
datum  surface  Is  calculated  by  the  fellswing  fortula  as  presented 
In  (6]: 

dB-a  (1  -  e2)  (A  (<p2  -  9^)  -  |_  (sin  2q>2  -  sin  2^) 

♦  £  (sin  4<p2  -  sin  ^)  -  (sin  6q>2  -  sin  69  ) 

♦  |  (sin  8q>2  -  sin  89^)-  (•*»  109g  -  Bin 

(83)  ♦  . ] 

The  parameters  In  the  above  equation  have  the  following  definitions: 
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a  »  ellipsoid  major  seal -diameter 
e  =  ellipsoid  eccentricity 


The  distance  between  two  points  on  a  parallel  on  the  datum 
surface  is  calculated  by  the  following  formula; 


(8*0  dp  -  N  cos  9  (12  -  XL) 

where  N  is  the  radius  of  curvature  in  the  prime  vertical. 

Utilizing  equations  (82),  (83),  and  (84),  the  linear 
distortions  along  meridians  and  parallels  at  a  particular  point  are 
calculated  by  the  following  ratios; 


(85) 


Meridian  Distortior. 


(86) 


Parallel  Distortion 


Equations  (85)  and  (86)  represent  ratios  of  projection  plane  distances 
to  datum  surface  distances.  By  making  the  distances  that  are  compared 
small,  the  above  equations  offer  a  good  approximation  of  tl. 
distortions  at  a  particular  point. 

5.2  Angular  Distortions 

The  meridians  and  parallels  of  a  sphere  or  an  ellipsoid 
of  revolution  intersect  at  right  angles  on  the  datum  surface.  When 
these  angles  are  projected  onto  the  projection  plane,  they  may  not 
intersect  at  right  angles  and  hence,  distortions  exist.  The  author 
chooses  to  indicate  this  distortion  by  the  absolute  value  of  the 
number  of  degrees  that  the  projected  angle  differs  from  ninety 
degrees . 

An  angle  on  the  projection  plane  can  be  calculated  from 
a  knowledge  of  the  coordinates  of  three  points. 


Y* 


Figure  l6 
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Figure  1 6  depicts  a  projection  plane  on  which  n  meridian,  a 
segment  of  which  is  P^P^,  and  '  parallel,  ’  segment  of  which  is 
PgPg,  intersect  at  P^.  From  Figure  l6,  it  is  evident  that  by 
knowing  the  coordinates  of  P^P^,  and  P2  in  a  Plane  coordinate 
system,  the  angles  0^  and  may  be  calculated  by  elementary 
geometry.  In  the  example  illustrated  by  Figure  l6,  the  angular 
distortion,  as  previously  defined,  has  the  value  (0  ♦  ®g) . 

Naturally  the  points  P  and  Pg  could  be  located  in  other  quadrants 
depending  on  the  projection,  but  the  basic  principles  in  determining 
the  angular  distortion  remains  the  same. 

Usually  the  meridians  and  parallels  are  projected  onto  the 
projection  plane  as  some  type  of  curved  lines.  Thus,  the  segment  of 
meridian  P^P^  and  segment  of  para  lied  PgP^  as  illustrated  in 
Figure  16  are  actually  curved  lines.  In  order  to  obtain  exact 
angular  values,  the  points  P^  and  P^  would  have  to  be  taken 
infinitely  close  to  Pq.  However,  it  will  be  illustrated  in  Chapter 
6  that  good  approximations  can  be  obtained  for  angular  distortions 

by  using  finite  distances  for  P  P  and  P  P  . 

10  2  0 

5.3  Distortion  Calculations 

The  methods  discussed  for  distortion  calculation  in  Sections 
5.1  and  5.2  indicate  that  basically  the  same  quantities  are  required 
to  calculate  both  linear  and  angular  distortions  at  a  particular 
point.  In  the  calculation  of  both  types  of  distortions,  the  plane 
X’-Y’  coordinates  of  three  points  are  required.  Since  the  calculation 
of  plane  coordinates  is  a  rather  lengthy  procedure,  the  distortion 


calculation  methods  discussed  in  this  paper  are  especia  ny 
for  use  when  electronic  computers  are  available. 


suited 
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6.  EXAMPLES  GP  PERSPECTIVE  MAP  PROJECTION  GRIDS 


6.1  General  Format  Of  The  Projection  Grids 

This  paper  has  dealt  with  the  derivation  of  General 
Perspective  Projection  Mapping  Equations  and  with  the  empirical 
derivation  of  the  distortions  inherent  with  the  perspective  pro¬ 
jection.  In  this  chapter  six  examples  of  the  perspective  projection 
are  Illustrated  and  briefly  discussed.  These  examples  are  Included 
both  to  Illustrate  various  types  of  perspective  projections  and  to 
partially  verify  the  validity  of  the  equations  derived  In  prior 
chapters.  The  coordinates  for  all  the  examples  were  computed 
through  the  use  of  equations  (78)  and  (8l).  All  calculations  were 
carried  out  on  the  IBM  709**  computer.  The  Appendix  Includes  the 
computer  program  used  to  produce  one  of  the  examples  and  also 
Includes  a  brief  discussion  of  the  program. 

A  similar  format  Is  used  for  all  the  Included  examples.  The 
first  page  of  each  example  lists  all  the  projection  parameters. 

The  body  of  each  example  Is  a  numerical  grid  which  tabulates  the 
coordinates  and  distortions  of  the  Intersections  of  meridians  and 
parallels  spaced  ten  degrees  apart.  The  grid  Is  read  by  selecting 
the  longitude  of  a  particular  desired  meridian  from  the  top  line 
and  then  by  reading  down  the  column  under  that  longitude  until  the 
desired  latitude,  as  Indicated  by  the  column  on  the  left  of  the 
grid.  Is  reached.  For  each  Intersection  of  a  meridian  and  a 
parallel,  five  quantities,  as  listed  on  the  left  of  the  grid. 
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are  tabulated.  The  first  two  quantities  are  the  X'  and  Y* 
coordinates  of  the  Intersection  as  computed  by  equations  (78)  and 
(8l) .  The  next  two  quantities  are  the  meridian  distortion  and  the 
parallel  distortion  computed  at  the  intersection  point  by  equations 
(85)  and  (86)  respectively.  The  fifth  quantity  is  the  angular 
distortion  of  the  intersection  of  the  meridian  and  the  parallel, 
calculated  as  described  in  Section  5 .2. 

Depending  on  the  type  of  projection,  the  grid  illustrates 
either  one  half  of  the  datum  surface  visible  from  the  projection 
center  or  a  grid  90  degrees  in  longitude  and  170  degrees  In 
latitude.  Im  both  of  the  above  cases  the  first  column  of  the  body 
of  the  grid  la  the  central  meridian  of  the  projection.  The  remaining 
columns  describe  meridians  to  the  right  of  the  central  meridian. 

The  grid  to  the  left  of  the  central  meridian,  If  calculated,  would 
be  iymmetrlcal  to  that  calculated  to  the  right  of  the  central 
meridian.  The  central  latitude  Is  found  in  the  center  of  the  grid 
and  may  be  located  by  the  latitude  Indicators  to  the  left  of  the  grid. 
tn  Inspection  of  these  latitude  indicators  reveals  that,  as  would 
be  expected,  the  absolute  values  of  the  latitudes  increase  until  a 
geographic  pole  Is  reached  and  then  they  proceed  to  decrease.  It  is 
evident  that  as  a  particular  meridian  passes  through  a  pole,  the 
meridian  undergoes  a  100  degree  change  In  longitude  and  after  pass* 
lng  through  the  pole  the  meridian  is  actually  to  the  left  of  the 
central  meridian.  Being  that  only  the  grid  to  the  right  of  the 
central  meridian  Is  being  considered,  as  a  meridian  passes  over  a 
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pole  the  symmetric  meridian  to  the  right  of  the  central  meridian  is 
considered  rather  than  the  original  meridian.  The  '■bore  can  be 
clarified  by  the  following  specific  example  of  a  meridian  passing 
over  a  pole.  In  the  examples  illustrated  in  this  chapter,  the  90 
Degrees  West  Longitude  Meridian  is  the  central  meridian.  Thus,  as 
the  80  Degrees  West  Longitude  Meridian  passes  through  a  pole,  it 
assumes  a  longitude  of  100  degrees  East  and  is  now  located  to  the 
left  of  the  central  meridian.  To  complete  the  right  side  of  the 
grid,  the  meridian  symmetric  to  the  100  Degrees  East  Longitude 
Meridian  is  considered.  This  symmetric  meridian  is  the  80  Degrees 
East  Longitude  Meridian.  Therefore,  the  longitude  Indicators 
located  at  the  top  of  the  grids  are  numerically  correct  but  may 
indicate  East  or  West  Longitude.  It  should  be  kept  in  mind  when 
viewing  the  grid  that  the  actual  value  of  a  meridian  in  no  way  enters 
into  the  computation  of  the  plane  X'-Y'  coordinates  of  points  on 
that  meridian.  It  is  the  longitude  difference  between  a  particular 
point  and  the  central  meridian  that  determines  the  values  of  the 
coordinates  of  that  point. 

It  will  be  noted  that  the  distortions  listed  for  90  degrees 
of  latitude  are  numbers  consisting  of  n  series  of  ones.  These 
numbers  are  not  the  distortions  at  this  latitude  but  merely 
indicators  that  the  distortions  at  this  latitude  were  not  calculated. 
It  is  apparent  from  the  distortion  analysis  in  Chapter  5  that  the 
equations  used  to  compute  distortions  for  the  rest  of  the  grid  will 
not  function  properly  at  the  pole.  No  effort  is  made  in  this  paper 
to  establish  special  distortion  equations  for  the  poles. 
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Earlier  in  this  paper  it  was  stated  that  to  calculate 
exact  distortions  an  infinitely  small  distance  between  two  points 
on  the  datum  surface  and  the  corresponding  distance  on  the  projection 
plane  would  have  to  be  compared.  The  distance  actually  selected 
would  determine  the  accuracy  of  the  calculations.  Various  small 
distances  were  selected  and  tested  by  calculating  distortions  at 
points  where  distortions  were  known  quantities.  It  was  determined 
that  by  using  a  distance  of  .00001  degrees  of  arc,  or  approximately 
1  meter  on  the  surface  of  the  earth,  errors  were  apparently 
eliminated  in  the  values  printed  out  in  the  examples.  The  manner 
in  which  the  length  of  the  distance  affected  the  accuracy  of  the 
distortion  computations  varied  according  to  the  type  of  projection. 

6.2  Sphere  And  Ellipsoid  Parameters 
The  sphere  that  was  selected  to  represent  the  datum  surface 
in  the  examples  is  only  one  of  many  possible  spherical  approximations 
of  the  shape  of  the  earth.  A  sphere  radius  of  6,371,224  meters  was 
selected.  This  radius  produces  a  sphere  with  approximately  the 
same  area  and  volume  as  the  International  Ellipsoid. 

The  International  Ellipsoid  was  selected  to  represent  the 
ellipsoidal  datum  surface  approximating  the  shape  of  the  earth.  The 
International  Ellipsoid  is  defined  as  an  ellipsoid  with  a  semi -major 
axis  equal  to  6,378,388  meters  and  a  flattening  equal  to  . 


6.3  Example  1  -  The  Gnomonic  Projection 
The  first  example  illustrates  a  Gnomonic  Projection.  Since 
the  projection  center  is  located  at  the  center  of  the  sphere  in  a 
Gnomonic  Projection,  it  is  obvious  thal  a  complete  hemisphere  can 
not  be  projected  onto  a  plane.  The  computer  was  programmed  not  to 
print  out  coordinates  larger  than  the  absolute  value  of  99,999>999> 
and  this  is  the  reason  that  a  complete  square  grid  was  not  printed. 
It  can  be  seen  that  both  the  size  of  the  coordinates  and  the  size 
of  the  distortions  increase  rapidly  as  the  location  of  points 
progress  away  from  the  central  meridian  and  the  central  parallel. 
The  property  of  a  Gnomonic  Projection  that  great  circles  on  the 
datum  surface  are  projected  as  straight  lines  is  clearly 
illustrated  by  the  projection  of  the  equator.  It  is  seen  that  the 
equator,  the  parallel  with  0  degrees  of  latitude,  has  a  constant 
Y'  coordinate,  readily  identifying  it  as  a  straight  line. 

The  following  five  pages  illustrate  a  Gnomonic  Projection. 
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PERSPECTIVE  PAP  PROJECTIONS. 


EXAPPLE  1  -  GNOMOMC  PROJECTION 


PROJECTION  PARAMETERS  - 


CENTRAL 

LATITULE  = 

40  CEGREES 

CENTRAL 

LONGITLDE  = 

9C  CEGREES 

HEIGHT 

CF  PROJECT  ION 

CENTER  =  - 

6371224.000 

PETERS 

FOCAL  LENGTH  (SCALE 

FACTOR)  =  - 

6371224. OOC 

PETERS 

SPHERE 

RADIUS  *  637122*4.000  PETERS 
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1.2376 

1.3716 

1.5896 

1.9654 

2.6511 

PAR. 

1.8273 

1.8199 

2.2377 

2.8655 

3.7559 

ANG. 

1.6960 

6.1398 

1  3. 0316 

22.1139 

32.7324 

© 
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LONG  I T UOC . . . . .  AC  3C  20  1C  C 


LAT. 

4C 

X 

47 30372. 32 

598 1928. P4 

74  7  C9II 8. 40 

9331603.41 

1 1812500. 7C 

V 

1417860.59 

2219975.14 

3362577.79 

50  3  3 1 20 . 1  8 

7592929.09 

MfcK. 

1 .296? 

1.4987 

1.0419 

2.4742 

3. 7653 

PAK. 

1  .5  758 

1.9413 

2.5104 

3.4431 

5.0984 

ANG. 

7 . 440  3 

13.0150 

20.8027 

3c'  .  7934 

42.3941 

LAT. 

30 

12445913.44 

17167853.75 

X 

5652046.28 

7316509.42 

94  56  387.  15 

Y 

2147Gb.  67 

1021258. 1 1 

2230371 .23 

4178814.84 

7592929.09 

MEK. 

l  .4479 

l  .  7542 

2. 3090 

3.4436 

6.2229 

PAR. 

1.6990 

2. 1825 

3.0100 

4.5727 

0.G424 

ANG. 

13.6920 

2i< .  35 1  1 

29.1908 

4C.2634 

52.9955 

LAT. 

20 

X 

6719319.00 

8943015.50 

1207156b. 38 

17097576.91 

27232625.23 

Y 

-1178*10.91 

-439665.70 

757559.56 

2902708.73 

7592929.09 

MER. 

1 . 7300 

2.2261 

3. 1958 

5.5198 

13.2993 

PAR. 

1 .9321 

2.6211 

3. 9482 

6.9752 

16.4093 

ANG. 

2C.5317 

28.4251 

38.4688 

50.7024 

64. 5862 

LAT. 

10 

X 

8057257. 36 

11116142.05 

15950704.66 

25468232.98 

56212979.68 

Y 

-2925073.40 

-2391562.69 

-14 38952. 71 

606583. 12 

7592929. G9 

MER. 

2.2754 

3.1315 

5 . 0802 

11.1511 

51.5931 

PAR. 

2.3515 

3.4487 

5.9224 

13.4624 

62. 1357 

ANG. 

28.3060 

37.5549 

48.8464 

62.1979 

77.03  75 

LAT. 

0 

X 

9911064.98 

14405539.75 

22850885. 78 

47160290.52 

Y 

-5346091 .71 

-534 6G 9 1.71 

-5346091 .71 

-5346091 .71 

MER. 

3.3396 

5 . 1 C  0  4 

10.1119 

37.0958 

PAK. 

3.1594 

5.22L6 

1 1 .1595 

43.2910 

ANG. 

37.4537 

48.0699 

60.4798 

74.6602 

LAT. 

-10 

X 

12875537. 14 

20459851 . 38 

40?  72557.57 

Y 

-9214090.95 

-10784059.71 

-15210889.76 

MER. 

5.8104 

10.6082 

32.3048 

PAR. 

4.9908 

10. 1 31 1 

34.2720 

ANG. 

48.4801 

60 . 2 1 8  3 

73. 3143 

LAT. 

-20 

X 

18084397.27 

37014592.43 

Y 

-17058900.76 

-2565  3407.89 

MER. 

13.7282 

38. 1  342 

PAR. 

1 0 . 7  o  7  6 

34 . 556  5 

ANG. 

b  1  .  H  1  8  3 

73.9726 

LAT. 

-30 

X 

4C239340. 38 

Y 

-44935799. 15 

MER. 

7  3 . 3  8  c>  0 

PAR. 

54 .03  1 1 

ANG. 

77.4316 

60 


f 

I 
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6.4  Example  2  -  The  Stereographic  Projection 
The  primary  property  of  a  Stereogrnphic  Projection  is  that 
it  is  a  conformal  type  projection.  Conformal  projections  have 
two  properties  which  are  immediately  apparent  in  the  following 
example.  First,  angles  are  projected  from  the  datum  surface  to  the 
projection  plane  with  true  size.  Thus,  it  is  seen  in  the  example 
that  all  angular  distortions  are  equal  to  zero,  or,  in  other  words, 
all  the  projected  meridians  and  parallels  intersect  at  right  angles 
The  second  property  inherent  in  conformal  projections  is  that  at 
any  particular  point,  scale  distortion  is  the  same  in  any  direction 
This  point  is  illustrated  in  the  example  by  the  fact  that  the 
meridian  and  parallel  distortions  are  equal  to  each  other  at  all 
intersections . 

The  following  five  pages  illustrate  a  Stereographic  Pro¬ 
jection. 
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PERSPECTIVE  MAP  PROJECTIONS 


EXAMPLE  2  -  STEREOGRAPHIC  PROJECTION 


• 

' 

— . - 

~ 

- — 

PROJECTION  PARAMETERS  - 

— 

— 

— 

CENTRAL  LATITUDE  *  40  DEGREES- 

CENT RAL  L ON G I XUD E  _* _ 90.  DEGRE E S _ _ 

HEIGHT  OF  PROJECTION  CENTER  »  - 12742448.000  METERS 


FOCAL  LENGTH  (SCALE  FACTOR)  *  -12742440.000  METERS 
$ PH ER E  RA D I  US  *  _ 637 1224.000  ME T E RS _ _ 


LONGITUDE .  90 

80 

70 

60 

?P 

LAT. 

50 

X 

V 

HER. 

-0.00 

12742448.00 

2.0000 

1411737.18 

12568441.75 

1.9851 

2720590.64 

12066612.20 

1.9423 

3841895.50 

11292150.08 

1.8762 

4721009. 34 
10321642. 19 
1.7934 

PAR. 

2.0000 

1.9851 

1.9423 

1.8762 

1.7934 

ANG. 

LAT. 

0.0000 

60 

0.0000 

0.0000 

0.0000 

o.oooo 

X 

-0.00 

938009.54 

1 B  2084 1 . 42 

2600577.24 

3241889. 38 

Y 

10692183.42 

10586682.43 

10279432.18 

9796365.59 

9175267.35 

HER. 

1.7041 

1.6957 

1.6712 

1.6327 

1.5832 

PAR. 

1.7041 

1.6957 

1.6712 

1.6327 

1.5832, 

ANG. 

LAT. 

0.0000 

70 

0.0000 

0.0000 

0.0000 

0.0000 

X 

-0.00 

562250.11 

1097776.10 

1582348.69 

1996255.98 

Y 

8922358.14 

8864353.79 

8694107.61 

8422399.70 

8065593.86 

MER. 

1.4903 

1.4859 

1.4729 

1.4523 

1.4252 

PAR. 

1.4903 

1.4859 

1.4729 

1.4523 

1.4252 

ANG. 

LAT. 

0.0000 

80 

0.0000 

0.0000 

0.0000 

O.OOOO 

X 

-0.00 

255809.90 

501842.04 

728907.43 

928925.75 

Y 

7356855.78 

7332571.55 

7260840.36 

7144931.85 

6989994.35 

HER. 

1.3333 

1.3315 

1.3262 

1.3177 

1.3062 

PAR. 

1.3333 

1.3315 

1.3262 

1.3177 

1.3062 

ANG. 

LAT. 

0.0000 

90 

0.0000 

0.0000 

0.0000 

0.0000 

X 

-0.00 

-0.00 

-0.00 

-0.00 

-0.00 

Y 

5941901.09 

5941901.09 

5941901.09 

5941901.09 

5941901.09 

HER. 

1.1111 

1.1111 

1.1111 

1.1111 

1.1111 

PAR. 

1.1111 

l.llll 

1.1111 

1.1111 

1.1111 

ANG. 

LAT. 

1.1111 

80 

1.1111 

1.1111 

1.1111 

1.1111 

X 

0.00 

217815.59 

4  30477 . 4  ~ 

632843.48 

819804.60 

Y 

4637871.78 

4655434.25 

4707826.  0«. 

479414B.25 

4912864.40 

HER. 

1.1325 

1.1338 

1.1376 

1.1440 

1.1528 

PAR. 

1.1325 

1.1338 

1.1376 

1.1440 

1.1528 

ANG. 

LAT. 

0.0000 

70 

0.0000 

0.0000 

0.0000 

0.0000 

X 

0.00 

406428.93 

805622.82 

1190157.15 

1552244.72 

Y 

3414328.65 

3444483.50 

3534796.53 

3684772.75 

389J451. 41 

HER. 

1.0718 

1.0741 

1.0809 

1.0923 

1.1082 

PAR. 

1.0718 

1.0741 

1.0809 

1.0923 

1.1082 

ANG. 

LAT. 

0.0000 

60 

0.0000 

0.0000 

0.0000 

0.0000 

X 

0.00 

572090.86 

1136958.30 

1686957.30 

2213600.73 

Y 

2246837.38 

2285770.51 

2402780.36 

2598445.63 

2873548.88 

HER. 

1.0311 

1.0342 

1.0435 

1.0591 

1.0810 

PAR. 

1.0311 

1.0342 

1.0435 

1.0591 

1.0810 

ANG. 

LAT. 

0.0000 

50 

0.0000 

0.0000 

0.0000 

0.0000 

X 

0.00 

719303.37 

1432848.83 

21 3428371 T 

”2816031750 

Y 

1114819.75 

1159488.58 

1294152.62 

1520744.24 

1842338.43 

HER. 

PAR. 

1.0077 

1.0115 

1.0230 

1.0423 

1.0697 

1.0077 

1.0115 

1.0230 

1.0423 

1.0697 

ANG. 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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LONGI  T  IJDE-  . _ 90 

80 

70 

6G 

50 

L.AT. 

AO 

X 

o.oo 

851309.16 

1699347.24 

2540174.00 

3368443.47 

Y 

0.00 

47874.76 

192605.37 

437505.40 

788066.07 

MER. 

1.0000 

1.0045 

1.OL80 

1.0439 

1.0737 

PAR. 

1.0000 

1 .0045 

1.0180 

1.0409 

1.0737 

ANG. 

LAT. 

0.0000 

30 

0.0000 

0.0000 

0.0000 

0.0000 

X 

0.00 

970389.81 

1940709.20 

2910260.49 

3876994.93 

Y 

-1114019.75 

-1065938.18 

-917792.49 

-665835.20 

-302348.17 

HER. 

1.0077 

1.0128 

1.0284 

1.0549 

1.0931 

PAR. 

1.0077 

1.0128 

1.0284 

1.0549 

1.0931 

ANG. 

LAT. 

0.0000 

20 

0.0000 

0.0000 

0.0000 

0.0000 

X 

0.00 

1078031.63 

2i59671.B8 

3248060.38 

4345290.57 

Y 

-2246837. 38 

-2198952. 12 

-2053495.86 

-1804966.80 

-1443860.14 

MER. 

1.0311 

1.0369 

1.0547 

1.0850 

1.1291 

PAR. 

1.0311 

1.0369 

1.0547 

1.0850 

1.1291 

ANG. 

LAT. 

0.0000 

10 

0.0000 

0.0000 

0.0000 

0.0000 

X 

0.00 

1174985.94 

2357536.27 

3555010.75 

4774264.56 

V 

-3414328.65 

-3369351.81 

-3232450.04 

-2997557.31 

-2654042.94 

MER. 

1.0718 

1.0784 

1.0986 

1.1332 

1.1838 

PAR. 

1.0718 

1.0784 

1.0986 

1.1332 

1.1838 

ANG. 

LAT. 

0.0000 

0 

0.0000 

0.0000 

0.0000 

0.0000 

X 

0.00 

1261225.87 

2534048.46 

3830209.56 

5161686.00 

Y 

-4637871.78 

-4597710.23 

-4475242.20 

-4264328.57 

-3954080.80 

MER. 

1.J325 

1.1400 

1.1629 

1.2023 

1.2604 

PAR. 

1.1325 

l . 1 400 

1.1629 

1.2023 

1.2604 

ANG. 

LAT. 

0.0000 

-10 

0.0000 

0.0000 

0.0000 

0.0000 

X 

0.00 

1335776.11 

2687026.06 

4069770.02 

5501130.51 

V 

-5941901.09 

-5908527.26 

-5806596.92 

-5630483.71 

-5370108.74 

MER. 

1.2174 

1.2260 

1.2521 

1.2973 

1 . 3640 

PAR. 

1.2174 

1.2260 

1.2521 

1.2973 

1.3640 

ANG. 

LAT. 

0.0000 

-20 

0.0000 

0.0000 

0.0000 

0.0000 

X 

0.00 

1396354.14 

2811600.79 

4265579.63 

5780112.95 

Y 

-7356855.78 

-7332360.25 

-7257449.91 

-7127679.27 

-6935021.58 

MER. 

1.3333 

1.3431 

1.3731 

1.4249 

1-5020 

PAR. 

1.3333 

1.3431 

1.3731 

1.4249 

1.5020 

ANG. 

LAT. 

0.0000 

-30 

0.0000 

0.0000 

0.0000 

0.0000 

X 

0.00 

1438694.61 

2898814.36 

4403054.66 

5976810.99 

Y 

-8922358.14 

-8908965.95 

-8867974.21 

-8796830.89 

-8690902.98 

MER. 

1.4903 

1.5016 

1.5361 

1.5960 

1.6852 

PAR. 

1.4903 

1.5016 

1.5361 

1.5960 

1.6852 

ANG. 

LAT. 

0.0000 

-40 

0.0000 

0.0000 

0.0000 

0.0000 

X 

0.00 

1455293.82 

2933037.90 

4457089.51 

605431 1.63 

Y 

-10692183.42 

-10692183.42 

-10692183.42  - 

-10692183.42 

-10692183.42 

MER. 

1.7041 

1.7171 

1.7571 

1.8264 

1.9298 

PAR. 

1.7041 

1.7171 

1.7571 

1.8264 

1.9298 

ANG. 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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LONG  1 1  UDf .  AO  30  20 _ 10 _ .£> 


LAT.  50 


X 

5335866.89 

5691969.65 

5813275.12 

5733356.59 

548925  1.42 

y 

9237041.51 

8112659.28 

7007  798. 45 

5964757.89 

5010423.49 

PLft. 

1 . 7008 

1.6049 

1 .5106 

1  .4216 

1.3401 

PAR. 

1.7008 

1.6049 

1.5106 

1.4216 

1. 3401 

ANG. 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

1  AT  . 

GO 

X 

1724348. 1 l 

4041756.60 

41 99426. 25 

4210577.24 

40^2853.57 

y 

8459536. 16 

7692281 . 59 

6911988.16 

6150125.56 

5430512.29 

HER. 

1 .5262 

1 .4650 

1 . 4029 

1.3421 

1.2848 

PAR. 

1.5262 

1 ,4650 

1.4029 

1.3421 

1.2848 

ANG. 

LAI. 

70 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

X 

2325529.35 

2562277.20 

2704245.76 

2753059.79 

2717027.44 

Y 

7643643.52 

7177963.50 

6 689545.93 

6197558.26 

5718499.97 

Mfc  R . 

1.3931 

1 . 3578 

1.3206 

1.2833 

1.2469 

PAR. 

1.3931 

1.3578 

1 . 3206 

1.2833 

1.2469 

ANG. 

1  AT 

80 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

X 

1095321.43 

1223264.16 

1 >09749.03 

t  353534.76 

1354974.  14 

Y 

6802652.21 

6590527.58 

6361 748.16 

6124494.27 

5886622.70 

HER. 

1.2924 

1.2767 

1.2598 

1.2423 

1.2247 

PAR. 

1.2924 

1.2767 

1.2598 

1.2423 

1.2247 

ANG. 

LAI. 

90 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

X 

-0.00 

-0.00 

-0.00 

-0.00 

-0.00 

Y 

5941901.09 

5941901 .09 

5941901.09 

5941901.09 

5941901.09 

HER. 

1.1111 

l.  1111 

l.llil 

l.llil 

l.llil 

PAR. 

1.1111 

l.llil 

l.llil 

l.llil 

l.llil 

ANG. 

LAT. 

00 

1 .1  111 

1.1111 

1.1111 

l.llil 

1.1111 

X 

986326. 36 

l 127519.46 

1238747.48 

1315777.39 

1  354974.  14 

Y 

5061746.98 

5237812.36 

54  37263,8  1 

5655386.65 

5886622. 70 

HER. 

1.1638 

1.1768 

1.1915 

1.2076 

1.2247 

PAR. 

1.1638 

l  .  1  768 

l.  I'M  4 

1 .2076 

1.2247 

ANG. 

LAT. 

70 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

X 

1083598. 1  2 

21 76352.68 

2418082.29 

2601948.39 

2717027.44 

Y 

4159200.46 

4  4  794  2  7.20 

4850210.15 

5265868.01 

5718499.97 

MLR. 

1.1284 

l.  152  7 

1.1809 

1.2125 

l • 2469 

PAR. 

1. 1284 

1.1527 

t . 1809 

1 .2125 

1.2469 

ANG. 
LAT  . 

60 

0.0000 

0  0000 

0.0000 

0.0000 

0.0000 

X 

2707147. 1 7 

3156217.94 

.354  7486.24 

3865513.96 

4092853.57 

Y 

3228781.84 

J66429 1 . 02 

4179028.18 

4  769072.43 

5430512.29 

HER. 

1.1093 

l  .1440 

1.1851 

1.2321 

1.2848 

PAR. 

1.1093 

1 . 1 440 

1.1851 

1.2321 

1.2848 

ANG. 

LAT. 

50 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

X 

346861 5.63 

40  79^(13.31 

46  34  305.  76 

5111988.83 

5488251.4 2 

Y 

2262892. 74 

2  786795.67 

3418128.94 

4159517.61 

5010423.49 

HER. 

1.  1056 

1.1503 

1.2042 

1.2675 

l. 3401 

PAR. 

1 . 1056 

1.1503 

1  .  2042 

l  .<.675 

1. 3401 

ANG. 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

65 


LONGITUDE.....  40  30  _ 20 _ IQ _ 0 


WSJj_ 

X 

V 

4176531,76 

1251660.20 

1*1171 

4953461.50 

1838296.93 

1.J719 

5683567.09 

2558086.76 

6344883.24 
3422194. 2B 

6907336.45 
4439950. ?9 

HER, 

1,1392 

1.3201 

1-4153 

PAR. 

ANG. 

LAT, 

30 

1.1171 

0.0000 

1.1719 

0.0000 

1.2392 

0.0000 

1.3201 

0.0000 

1.4153 

0.0000 

X 

Y 

HER. 

4836564.26 

183729.37 

1.144? 

5781157.52 

806949.55 

1.2098 

6697545.62 

1585340.49 

1.2917 

7564058.05 

2539961.51 

JL.J922 

8351245.20 

3693555.04' 

1.5136 

PAR. 

ANG. 

LAT. 

20 

1.1443 

0.0000 

1.2098 

0.0000 

1.2917 

0.0000 

1.3922 

0.0000 

1.5136 

0.0000 

_  X 

V 

HER. 

5451593.77 

-956169.27 

1.1887 

65  .i4106.26 
-322711.33 
1.2660 

7674962.10 

481647.59 

1.3642 

8768342.18 

1488669.71 

1.4872 

9815977.84 
2736865. 18 
L,5396 

PAR. 

ANG. 

LAT. 

1.1687 

0.0000 

»(> 

1.2660 

0.0000 

1.3642 

0.0000 

1.4872 

0.0000 

1.6396 

0.0000 

X 

Y 

HER. 

6021129.02 

-2185885.79 

1,2527 

7299483.63 

-1570434.47 

1.3433 

8609607.63 

-776694.10 

1.4602 

9945288.24 

236869.36 

J.6095 

11288816.35 

1524829.01 

1.7992 

PAR. 

ANG. 

LAT. 

5 

1.2527 

0.0000 

1.3433 

0.0000 

1.4602 

0.0000 

1.6095 

0.0000 

1.7992 

0.0000 

X 

6540643.34 

7979108.00 

9488081.93 

11075565.28 

12742448.00 

Y 

-3527780.01 

-2961155.49 

-2219789.49 

-1255313.41 

0.00 

HER. 

1.3401 

1.4461 

1 . 5848 

1.7652 

2.0000 

PAR. 

1.3401 

1.4461 

1.5848 

1.7652 

2.0000 

ANG. 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

LAT. 

-jo 

X 
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6.5  Example  3  -  The  Orthographic  Projection 
Rather  than  illustrating  an  Orthographic  Projection  with  the 
projection  center  located  at  an  arbitrary  latitude  of  Uo  degrees  as 
was  done  in  Examples  1  and  2,  a  Polar  Orthographic  is  illustrated 
in  Example  3.  A  Polar  Orthographic  is  a  perspective  projection  in 
which  the  projection  axis  is  coincident  with  the  axis  of  rotation 
of  the  datum  body,  and  the  projection  center  is  located  at  an  infinite 
distance  along  this  axis.  In  this  type  of  projection  the  most 
prominent  properties  are  that  all  parallels  of  latitude  are  projected 
in  true  size  as  concentric  circles  with  centers  at  the  pole  "nd  nil 
meridians  are  projected  as  straight  lines  radiating  from  the  pole. 

These  properties  are  illustrated  in  the  following  example  by  the 
parallel  distortions  all  being  equal  to  unity,  or,  in  other  words, 
the  distances  along  the  parallels  measured  on  the  projection  plane 
are  equal  to  the  distances  along  the  parallels  measured  on  the 
datum  surface.  Also,  it  is  noted  that  all  the  angular  distortions 
are  equal  to  zero.  This  is  due  to  the  fact  that  the  parallels  are 
all  concentric  circles  and  the  meridians  are  all  straight  lines 
radiating  from  the  center  of  the  concentric  circles.  Thus,  all 
intersections  of  meridians  and  parallels  form  right  angles. 

The  following  five  pages  illustrate  an  Orthographic  Projection. 
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PERSPECTIVE  MAP  PkOJCCIIONS 
EXAMPLE  3  -  OR  f  HUGP  APH i C  PROJECTION 
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_  .  C..OC 

HER. 

0.6428 

0.6428 

0.6428 

0.6428 

6.6428 

PAR. 

l.OlOO 

1.0000 

l.OOCJ 

1  .COCO 

1.0000 

ANG. 

0  .O  JDO 

O.OOCC 

0.0000 

0.0000 

C.CCOf 

LAT. 

30 

X 

4226758.87 

4770418.00 

5184887. 32 

5433816.46 

5517641.84 

Y 

-3546671.81 

-2  75b8  20.92 

-  1887144.65 

-958128.45 

0.00 

HER. 

0.5000 

G.5CCC; 

0.50CG 

G.500C 

o.5rcc 

PAR. 

1.0000 

1 .0000 

1.0000 

l.COCO 

l.OOOC 

ANG. 

P.C'jOG 

C .COOC 

o.oocc 

0  •  C  0 1  ■  0 

o.oocc 

LAT. 

20 

X 

4586302.09 

5184887. 32 

5625932. 37 

5896C  36.31 

5986992.  18 

Y 

-3848364. 39 

-299 3496.09 

-2047671.9? 

-103963C.2B 

c.oc 

HER. 

0.3420 

0.3420 

0. 3420 

0.3420 

C. 342C 

PAR. 

1 .0'  Ou 

1 .coon 

1 .0000 

l.CCOO 

l.OOCC 

ANG. 

0  •  0  C>  0  0 

O.OOCC 

0.0000 

G.GCCO 

0.0000 

LAT  . 

10 

X 

4006492.84 

54  338  16.4*. 

'*896036.31 

61  791C8.C9 

6274430. 79 

Y 

-4033126.37 

-313 1/ 15.40 

-2145981.7? 

-1089543.47 

0  •  •  1 0 

HER. 

0.1736 

0.1736 

0.1736 

0.1736 

c. 1736 

PAR. 

l.OOCO 

l .0002 

l.CCOO 

1  .OCCG 

l.GC  0 

ANG. 

y*  #0'.  00 

0.0*00 

0,0000 

G.COOO 

0.000 o 
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6.6  Example  h  -  A  General  Projection  Of  A  Spherical  Dntum 

In  Examples  1,  2,  and  3,  the  parameter  h,  thp  height  of  the 
projection  center  above  the  datum  shpere,  was  set  equal  to  a 
particular  value  in  ielationship  to  the  datum  sphere  in  order  that 
the  resulting  projection  would  exhibit  certain  desired  properties. 

In  Example  U,  a  completely  arbitrary  value  o  approximately  700 
miles  has  been  selected  for  the  parameter  h.  Thus,  ideally.  Example 
k  illustrates  the  type  of  projection  that  would  be  obtained  by 
taking  a  vertical  photograph  of  the  datum  sphere  from  a  height  of 
approximately  700  miles. 

In  the  following  example  the  computer  was  programmed  to  prirt 

out  the  coordinates  of  only  those  meridian-parallel  intersections 

that  are  visible  from  the  projection  center.  By  visible,  the 

author  means  that  the  intersection  point  is  not  beyond  the  horizon 

as  viewed  from  the  projection  center.  The  relationship  between 

latitude  and  longitude  values  on  the  horizon  is  calculated  by 

setting  LP  equal  to  0  in  equation  (23),  which  results  in  the 
I 

following  equation: 

(87)  cos  9  ■  cos  t  cos  $  cos  \  +  sin  |  sin  $ 

c  c 

It  is  interesting  to  note  that  the  shape  of  the  projected  grid  in 
the  following  example  actually  resembles  the  shape  of  the  grid  th”t 
would  be  visible  from  the  projection  center. 

The  following  two  pages  illustrate  a  general  projection 
of  a  spherical  datum. 
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PERSPECTIVE  MAP  PROJECTIONS 


EXAMPLE  A  -  GENERAL  PROJECTION 


PROJECTION  PARAMETERS  - 

“CENTRAL  LATITUDE  *  40  DEGREES- 

_ CENTRAL  LONGITUDE  *  90  DEGREES _ 

HEIGHT  OF  PROJECTION  CENTER  =  1126542.900  METERS 

FOCAL  LEN  G T El  '(SCALE  FACTOR I~  Tf26542 7900  METERS 
SPHERE  RAO  I  US  *  6371224. GOO  METERS 


LONGITUDE.....  90  80  _  -70 - 60 - 54 


1  AT 

70 

X 

0.00 

212555.86 

4C3501 . 10 

Y 

1812374.71 

1801410.01 

1770424.38 

HER. 

0.0351 

0.1495 

0.2695 

PAR 

0.6689 

0.5421 

0.4696 

ANG. 

0.0000 

79.6911 

89.7720 

LAT  . 

60 

X 

0.00 

402607.95 

740325.02 

976410.19 

1108127.44 

Y 

1624884.93 

1608607.03 

1564559.0/ 

1503979.25 

1438497.17 

HER. 

0.3328 

0.3861 

0.4715 

0.5142 

0.5060 

PAR. 

ANG. 

LAT. 

60 

0.7457 

0.0000 

0.6845 

35.2465 

0.5306 

56.0944 

0.3475 

70.6413 

0. 1920 
89.7929 

X 

0.00 

630323.72 

1  117097.50 

1403461.63 

1515120.17 

Y 

1018814.33 

1016056.92 

1008965.24 

1000010.81 

991240.37 

HER . 

0.7623 

0.7612 

0. 7354 

0.6673 

0.5740 

PAR 

0.9209 

0.8069 

0.5424 

0.2668 

0.0608 

ANG. 

0.0000 

15.7277 

26.8938 

33.7398 

43.8255 

LAT. 

<.0 

X 

6.00 

806833.61 

1390891. 7e 

1689227.46 

1765997.41 

Y 

0.00 

45373.60 

157644. 7B 

290943.11 

413164.91 

HER. 

1.0000 

0.9517 

0.8301 

0.6833 

0.5482 

PAR. 

1.0000 

0.8529 

0.5357 

0.2539 

0.1310 

ANG. 

LAT. 

30 

0.0000 

0.7103 

5.9601 

24.3073 

74.7379 

X 

0.00 

838315.20 

1438161.72 

1736647.90 

Y 

-IC;8814.33 

-920859.00 

-680129.73 

-397325.71 

HER. 

0.7623 

0.7302 

C. 6481 

0.5463 

PAR. 

0.9209 

0.7987 

0. 5438 

0.3357 

ANG. 

LAT. 

20 

0.0000 

17.1789 

37.3981 

65.2857 

X 

0.00 

741046.79 

129061 1.26 

Y 

-1624884.83 

-1511575.70 

-1227160.90 

HER. 

0.3328 

0.3572 

0.3854 

PAR. 

0.7457 

0.6722 

0.5137 

ANG. 

LAT. 

10 

0.0000 

37  6277 

66.9396 

1 

0.00 

5  i  7-122.64 

i 

-1812374.71 

-l 71  296.95 

HER. 

0.0151 

0.1408 

PAR. 

0.5689 

G. 4304 

ANG. 

o.oooc 

95.9367 

6.7  Example  5  -  A  general  Projection  Of  An  Ellipsoidal  Datum 

Example  5  illustrates  a  projection  analogous  to  Example  4 
with  the  exception  that  an  ellipsoidal  datum  surface  is  utilized 
rather  than  a  spherical  datum  surface.  The  same  value  of  h, 
approximately  700  miles,  is  also  used  in  Example  5.  Therefore, 
Example  5  illustrates  to  a  better  approximation  than  Example  4 
the  projection  of  the  actual  meridians  and  parallels  on  the  surface 
oi  the  earth. 

Equation  (87)  was  again  utilized  in  obtaining  those 
intersections  of  meridians  and  parallels  that  are  visible  from  the 
projection  center.  While  equation  (87)  is  an  exact  formula  when 
dealing  with  a  spherical  datum  surface,  this  equation  is  only  an 
approximation  when  dealing  with  ar  ellipsoidal  datum  surface.  The 
author  wishes  to  emphasize,  however,  that  the  calculated  values  are 
exact  and  the  approximation  only  involves  the  determination  as  to 
which  intersections  are  visible. 

The  followini  two  pages  illustrate  a  general  projection  of 
an  ellipsoidal  datum. 
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PERSPECTIVE  MAP  PROJECTIONS 
EXAMPLE  5  -  GENERAL  PROJECTION 


PROJECTION  PARAMETERS  - 

~~  CENTRAL  LATITUDE  *"  40  DEGREES 

_ CENTRAL  LONGITUDE  =  90  DEGREES _ _ _ 

HEIGHT  OF  PROJECTION  CENTER  =  1126542.900  METERS 

FOCAL  L ENG TH  (  SCALE  FACTOR)  =  l  1^3427900  METERS 

_ ELLIPSOID  MAJOR  SEMI  DIAMETER  =  6378388.000  METERS 

ELLIPSOID  FLATTENING  =  1/297.0 


LONGITUOE. 

....  90 

80 

_  .7  Cl _ 

_  6C 

LAT. 

70  _ 

X 

O.Ou 

213267.66 

606766. 70 

Y 

1812769.58 

1801766.18 

1770668.50 

MER. 

0.0352 

0. 15C7 

0,2716 

PAR. 

0.5685 

0.5616 

0.6690 

ANG. 

0.0000 

79.7639 

89.8051 

LAT. 

60 

X 

0.00 

603959.06 

742663.99 

979168.62 

Y 

1626973.58 

1608639.02 

1 566650.96 

1503710. 6C 

HER. 

0.3366 

0.3881 

0.6760 

__0*J>167 

PAR. 

0.7655 

0.6862 

0.530C 

0.3466 

ANG. 

0.0000 

35.2651 

56. 1 362 

70.6878 

LAT. 

50 

X 

0.00 

632220.67 

1120121.61 

1 406725.49 

Y 

1018253.01 

1015508.62 

1008652.65 

999549.64 

MER. 

0.7637 

0.7626 

0.7366 

0.6600 

PAR. 

0.9209 

0.8067 

0.5616 

0.2658 

ANG. 

0.0000 

15.7602 

26.9061 

33.7372 

LAT. 

60 

X 

0.00 

800767.25 

1393808.42 

1692169.17 

V 

-0.00 

65682.36 

157975.36 

291449.78 

MER. 

1.0000 

0.9516 

0.8297 

C.6827 

PAR. 

1.0000 

0.8525 

6.5349 

1.2531 

ANG. 

O.OCOO 

0.7120 

5.975? 

24.39C9 

LAT. 

30 

X 

0.00 

860063.66 

1440743.42 

1739222.20 

Y 

- 

1016690.79 

-918639.66 

-677765.04 

-394953.84 

MER. 

0.7627 

C . 7  306 

0.6483 

6.5463 

PAR. 

0.9211 

0.7985 

0.5431 

C.334B 

ANG. 

0.0000 

17.1835 

37.4138 

65.3518 

LAT. 

20 

X 

0.00 

762761.38 

1293251.00 

V 

- 

1621626.01 

-1508007.55 

-1223424.93 

MER. 

0.3339 

0.3586 

0. 3865 

PAR. 

0.7663 

0.6725 

0.5133 

ANG. 

O.OCOO 

37.5217 

66.9340 

LAT. 

10 

X 

0.00 

599603.56 

V 

- 

1809173.66 

-1  710917.70 

MER. 

0.0359 

0.1618 

PAR. 

0. 5700 

0.5311 

ANG. 

0.0000 

85.6795 

_.50 


1110890. 16 
1438096.83 
0*50,82 
0.  1912 
89.9032 


1518C77. 95 
990836.13 
766 
O.CbCC 
63.8693 

176869  7.1.6 
613769.83 
C. 5675 
0.13C7 
76.6825 
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6.8  Example  6  -  An  Orthographic  Projection  Of  An  Ellipsoidal  Datum 

Example  6  is  included  as  a  final  illustration  of  the  utility 
of  the  mapping  equations  derived  previously  in  this  p^per.  The 
Orthogr.  phic  Projection  is  a  rather  commonly  used  projection,  hut 
to  the  best  of  the  author's  knowledge  a  spherical  d^tum  surface  is 
always  used  for  this  projection.  Example  6  illustrates  a  Pol«r 
Orthographic  Projection  using  an  ellipsoidal  datum  surface.  This 
type  of  projection  more  closely  approximates  the  projection  of  the 
actual  meridians  and  parallels  of  the  earth  than  a  projection  using 
a  spherical  datum  surface.  Therefore,  in  this  respect,  it  appears 
to  he  a  superior  type  projection. 

It  can  he  seen  that  the  properties  noted  in  Example  3  are 
alLO  evident  in  Example  6.  As  in  Example  3,  the  linear  distortions 
along  the  parallels  and  the  angular  distortions  in  the  intersections 
of  meridians  and  parallel  are  both  equal  to  zero. 

A  d.‘ sadvantage  in  using  an  ellipsoidal  datum  surface  rather 
than  a  spherical  datum  surface  is  obviously  the  more  complicated 
mapping  equations  required  for  the  former.  However,  as  the  use  of 
electronic  computers  increases,  this  disadvantage  correspondingly 
decreases. 

The  following  five  pages  illustrate  an  Orthographic 
Projection  of  an  ellipsoidal  datum. 
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PER  SPEC T 1 Vt  MAP  PROJECTIONS 


EXAMPLE  6  -  ORTHOGRAPHIC  PROJECTION 


PROJECTION  PARAMETERS  -  _ _ _ _ _ 

CENTRAL  LATITUDE  =  00  DEGREES 

C  EN  TR AL  LONG  I TUUC~=  90~UEGR?ES 

_ HEIGHT  OF  PRO JEC T i ON  CENTER  «  1 NF INI TY _ 

FOCAL  LENGTH  =  INFINITY 

ELL  IP  SO  ID  MAJOR  SEMI  ~ D I A  MfcT¥R~=  6378388.000  METERS 

_ _ ELLIPSOID  F  LATTEN  I NC,  =  1/29  7.0__  _ _ 


LONG  1 T  UOE 

HO 

70 

6C 

50 

LAT  . 
X 

V 

MLR. 

PAR. 

ANG. 

0 

-O.CO 

6370300.00 

1  IC7595.45 
6281486.95 

2181537.18 
5993  724.14 

3189194. OC 
5523846. C4 

4C99948. 78 
4886128.68 

C.O  '00 

1 .OGCO 

0  .Ot-OG 

(  .  C  GC  C 

1  .COCO 
C.GOCC 

O.COCC 

l.COCf 

c.rocc 

(.  .OCf  0 
l.CCCO 

c.ooco 

0.09GC 

l.CCCO 

C.0000 

LAT  . 
X 

V 

10 

-0.00 
6202122. 72 

1C90079. 16 
6186683. 16 

2148612.51 

5903264.36 

3141661. 36 
5440477.87 

4038070.65 

4812385.20 

MER. 

PAR. 

0.1736 

1  .C.CO 

C. 1 736 

l.coo: 

0.1736 

l.OCOO 

C. 1 736 
1.0000 

0.1736 

l.OOCC 

ANC  . 

6.0000 

0.0  600 

O.CCOO 

O.COCC 

C.OOOC 

LAT 

X 

V 

20 

-0.00 

5996082.27 

1041208.76 

5904988.31 

2C5C780.92 

56  34  4  74.26 

2998C41.13 

5192759.57 

3854207.39 
4593265. 50 

HER. 

PAR. 

ANG. 

C  .3420 
1.0(00 
0.0-00 

G. 3429 

0. 3420 
l.OCOO 
0.0000 

C.342C 

0.3420 

1 .0003 
0.0009 

1 .0000 
C.COOG 

l.OOCC 

O.OOCG 

LAT  . 

X 

y 

30 

-0.00 

5528493.78 

960012.07 

5444503.53 

1890866.23 

51  95c84 . 8  1 

2764246.89 

4787816.06 

3553647.30 

4235C71.94 

HER. 

0.5CGO 

0.5CCC 

1  .CO 30 

o.rcco 

C.50C0 

1.0000 

0.0000 

C.5C0C 

0.5000 

PAR. 

ANG. 

l.OCOO 

O.CCOO 

l.CCCO 

c.cooe 

l.OCOO 

O.OOGO 

LAT. 

X 

Y 

40 

-G  .00 
4892928.82 

04964H. 1 7 
4818594.24 

1673460.22 
4597849. 1 l 

2446464.41 

42374CC.66 

3145114.02 

3748203.93 

MER. 

PAR. 

ANG. 

C .6428 

1 .0000 

0  .uOOO 

0.6428 
1.0000 
’ .  G  C  C  ^ 

0.6428 

t.coco 

C.COOo 

L .6428 
l.COCC 
C.OOCO 

0.6426 

l.OCOO 

O.OCGO 

LAT. 

50 

X 

-O.CO 

71 3357. 1 3 

1405039.26 

2954029.99 

2640610.06 

Y 

4108359.98 

4C45649. 32 

386031 3.65 

3557684.31 

314695ft. 52 

MER. 

0.7674 

C.7674 

0.7674 

G.  76  74 

0.7674 

PAR. 

l.OvCJ 

l  t  C  L'lJ 

l.OOCC 

1  .CCCC 

l.OC JC 

ANG. 

C.OOCO 

O.OCGO 

o.occo 

:  .core 

O.CCOO 

LAT. 

60 

X 

-0.00 

556199.15 

109352H. 66 

1590632.25 

2055162. JO 

Y 

3197264.49 

1  14 8690.86 

3004445.85 

2768912.27 

2449246. 70 

MCR. 

3.87C4 

C.O  704 

0.0704 

'  .0  704 

C.87C4 

PAR. 

1.0. C  ! 

1  .COCO 

l.O'CC 

l  .001*0 

1.09)9 

ANG. 

0.0  0() 

0.0006 

o.cccc 

C.COOO 

C.0301 

1  AT  . 

70 

X 

-O.oO 

379949.37 

748  334. 18 

1094'-  2C.62 

1406445. a: 

Y 

2188341.24 

2 l 54799. 78 

2056966.21 

l  994899.30 

1676136.84 

ME  R. 

r .9473 

G  .  94  7  C 

0.9470 

C  .94  70 

0.94  7C 

PAR. 

1  .  OgC 3 

l  .  C  3C  0 

1.0600 

l.CCCO 

1 .  O'*  L 

ANG. 

0.0-00 

>-.COCO 

O.OCOC 

f.  .CCCC 

C .00  cc 

LAI  . 

00 

X 

-  j.90 

l  ’*  2 '362  .  (.  1 

JUCC60.97 

555611.97 

714280.98 

Y 

1111221.94 

1 094 341 .95 

1044268.94 

962348.16 

851246. 93 

MER. 

6.99 42 

C .9942 

0 .9942 

V.-.994  2 

*.9942 

PAR. 

l  .coco 

i  3c  c 

1.3  30C 

l  .ooco 

1.0000 

ANG  . 

C  .  Ji.CO 

c . C  GcC 

O.OOCC 

c.ooco 

O.OCCl 

81- 


LONGITUDE 

80 

70 

60 

5C 

LAI  • 

9C 

X 

-3.00 

-C.C3 

-O.CC 

-O.CC 

“  C  •  00 

V 

3.00 

-C.t3 

-C.00 

-:.oc 

-c.c: 

HER. 

1.1111 

1.1111 

1.1111 

1.1111 

l.iiu 

PAR. 

1.1111 

1.1111 

1.1111 

1.1111 

l.iiu 

ANG. 

1.1111 

1.1111 

1.1111 

1.1111 

l.iiu 

LAT. 

8G 

X 

0.30 

192962.01 

380060.97 

555611.97 

716280.98 

V 

-1111223.96 

. -1096J61.95 

-1066208.96 

-962368.16 

-851266.93 

HER. 

C.9962 

0.996? 

0.9962 

0.9962 

C.9962 

PAR. 

1.0030 

I.U3C0 

l.COGC 

1.0000 

1.03G0 

ANG. 

C .0303 

c.occc 

O.OCOO 

C.OCCC 

fa.OGCfa 

LAI. 

70 

X 

3.33 

379969.37 

768)56.18 

1096C2C.62 

1606665.83 

V 

-2188361.26 

-2156799.98 

-2C56C86.21 

-1896899.30 

-1676136.86 

HER. 

f .9673 

0.9673 

3.9670 

W.967C 

0.967C 

PAR. 

1.GC39 

1  OOOC 

1.0000 

l.COCC 

1 .00 vc 

ANG. 

C.OGOO 

C.OOCC 

0.0000 

C.OOOC 

C.00 oC 

LAT. 

60 

X 

0.30 

555199.15 

1093528.86 

1598632.25 

2055 162. CO 

V 

-3197266.69 

-316R69C.86 

-3C06665.H5 

-2768912.27 

-2669266.  7** 

HER. 

3. 9706 

C.87'6 

0.8706 

C.87C6 

0.8736 

PAR. 

1.0CC3 

1.33L0 

l.COCC 

l.CCCC 

1.00 ;c 

ANG. 

C.OCCO 

O.GOCC 

0. )OCC 

c.cooc 

c.ooj: 

LAT. 

5C 

X 

C.00 

713357.13 

1605039.26 

2C56C29.99 

2660610. C  6 

V 

—6108059.98 

-6C6S669. 32 

-3860)13.65 

-3557686.31 

-3166956.52 

HER. 

0.76/6 

0.7676 

0.7676 

0.76/6 

G. 7676 

PAR. 

1.0300 

1.00C9 

l.CCCC 

l.CCCC 

l.OCCC. 

ANG. 

_ 3.0300 

C.onoc 

C.OCCC 

C.OCCC 

o.eocc 

LAT. 

60 

X 

3.03 

869660.17 

1673680.22 

2666666.61 

3165116.02 

V 

— 6R9292B. 8? 

-6818596.26 

-6597869.11 

-62376CI .66 

-3768200.93 

HER. 

0.6628 

C.6628 

C.662B 

C.6628 

* .6926 

PAR. 

1 .G  CG3 

1 .OGCC 

l.CCCC 

l.CCCC 

l.COO' 
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7.  SUMMARY  AND  CONCLUSIONS 


The  primary  goal  of  this  paper,  as  stated  In  the  Introduction, 
was  to  develop  general  mapping  equations  for  the  perspective 
projection  which  would  be  suitable  both  for  spherical  and 
ellipsoidal  datum  surfaces.  A  careful  method  of  proceeding  from 
the  known  to  the  unknown  was  utilized  in  developing  the  desired 
equations.  First,  mapping  equations  for  selected  perspective  pro¬ 
jections  of  a  spherical  datum  surface  were  derived  utilizing 
already  validated  methods.  Next,  mapping  equations  for  a  spherical 
datum  surface  were  derived  by  a  novel  photogrammetrlc  method. 

These  mapping  equations  were  proved  equal  to  thoae  derived  earlier. 
After  having  at  least  partially  proved  the  photogrnmmetilc  method  •:.< 
be  valid  by  Its  use  on  a  spherical  datum  surface,  this  method  was 
used  ko  derive  mapping  equations  for  an  ellipsoidal  d«tum  surface. 
Equations  (78)  and  (8l),  the  General  Perspective  Projection  Mapping 
Equations,  were  the  results  of  the  derivation.  Equations  (78)  and 
(8l)  were  proved  to  be  general  equations  encompassing  all  the 
special  cases  that  had  been  discussed  prior  to  their  derivation. 

It  Is  evident  from  the  derivation  of  these  equations  that  they  could 
be  made  to  be  even  more  general  by  providing  for  the  possibility  that 
the  arbitrary  point  P  that  Is  to  be  projected  is  not  on  the  datum 
surface  but  at  some  distanc.  above  or  below  it.  The  alteration 
to  equations  (78)  and  (8l)  to  take  into  consideration  this  distance 
would  be  accomplished  by  first  chanrMnrr  coordinates  of  point 


85 


P,  as  given  "by  equations  (62),  (63),  and  (64),  to  "ccount  for  the 
distance  and  then  by  using  these  new  coordinates  in  the  remainder  of 
the  derivation.  Finally,  ns  a  means  of  vizunlizing  certain 
perspective  projections,  an  empirical  method  was  developed  to 
compute  various  distortions  in  the  projections  and  examples  of 
several  projected  grids  were  illustrated. 

The  author  believes  that  equations  (78)  and  (8l)  derived 
in  this  thesis  are  nseful  equations  in  that  they  are  truly 
General  Perspective  Projection  Mapping  Equations  suitable  for  both 
spherical  and  ellipsoidal  datum  surfaces  that  can  be  used  regardless 
of  the  location  of  the  projection  center. 
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APPENDIX 


The  computer  programs  utilized  to  produce  all  examples 
included  in  this  thesis  were  basically  very  similar.  The  actual 
equations  used  to  calculate  the  desired  data  were  identical  in 
all  the  examples,  and  the  primary  differences  in  the  programs 
concerned  format.  Thus,  although  this  appendix  is  devoted  to  an 
explanation  of  the  computer  program  used  to  produce  Example  5>  it 
essentially  explains  the  computer  programs  used  for  all  the 
Included  examples.  The  IBM  709^  computer  of  the  Ohio  State 
University  was  employed  to  produce  the  examples;  and  the  pro¬ 
gramming  language  used  was  SCATRAN,  a  program  language  taught  and 
used  at  Ohio  State  University, 

A  complete  and  detailed  explanation  of  the  program  will  not 
be  given,  but  sufficient  explanation  will  be  given  to  allow  the 
reader  to  understand  the  general  principles  of  the  program.  A 
more  thorough  discussion  of  the  computer  program  is  not  being 
presented  due  to  the  authors  belief  that  the  program  is  not  an 
end  product  but  only  a  tool  that  was  used  to  illustrate  and  verify 
the  concepts  developed  in  this  thesis. 

In  succeeding  paragraphs,  first  the  major  terms  used  in  the 
program  will  be  identified,  and  then  the  functions  of  the  Source 
Language  Statements  will  be  explained.  The  program  itself  is 
included  at  the  end  of  this  appendix. 

The  following  is  a  lif  of  the  major  terms  th"t  were 
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components  of  the  program  used  to  produce  Example  5  of  this 
thesist 

X  -  the  abscissa  or  X  coordinate  of  a  particular  point  In  meters 

V  -  the  ordinate  or  Y  coordinate  of  a  particular  point  in  meters 

V  -  a  longitude  value  used  in  printing  out  the  grid  format 

V  -  a  latitude  value  used  in  printing  out  the  grid  format 

IMP  -  an  angular  distortion  in  the  intersection  o**  °  meridian  and 
a  parallel 

RDM  -  a  ratio  of  a  meridian  distance  on  the  projection  plane  to  a 
meridian  distance  on  a  datum  surface 
RDP  -  a  ratio  of  a  parallel  distance  on  the  projection  plane  to  a 
parallel  distance  on  a  datum  surface 
RAD  -  a  conversion  factor  to  convert  degrees  to  radians 
A  -  the  sphere  radius  or  ellipsoid  semi -major  axis  in  meters 
FLAT  -  the  ellipsoid  flattening 

EC  -  the  height  of  the  projection  center  above  the  datum  surface 
in  meters 

F  -  the  distance  between  the  projection  center  and  the  projection 
plane  in  meters 

BC  -  the  central  latitude  in  degrees 
BP  -  the  latitude  of  a  particular  point  in  degrees 
I £  -  the  central  longitude  in  degrees 
LP  -  the  longitude  of  a  particular  point  in  degrees 
DBQCH  -  a  small  number  of  degrees  used  to  locate  a  point  very  close 
to  a  primary  point.  In  the  following  definitions  the 
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primary  point,  which  is  the  intersection  of  a  particular 

meridian  and  a  particular  parallel,  will  be  identified 

as  P  .  Figure  l6  illustrates  P  and  its  relationship 
0  0 

to  P,  and  P  . 

1  2 

BPRB  -  the  latitude  of  P^  in  radians 
LPRL  -  the  longitude  of  P2  in  radians 

U  -  a  term  used  to  indicate  the  computation  of  an  Orthographic 

Projection,  When  an  Orthographic  Projection  is  being  calculated, 
U  is  set  equal  to  0  and  F  is  set  equal  to  1. 

E2  -  the  square  of  the  eccentricity  of  an  ellipsoid 
Al,  Bl,  Cl,  Dl,  El,  FI  -  terms  used  in  equation  (83) 

IIP  -  the  radius  of  curvature  in  the  prime  vertical  at  a  particular 
point 

NC  -  the  radius  of  curvature  in  the  prime  veritcal  along  the 
projection  axis 

NPB  -  the  radius  of  curvature  in  the  prime  vertical  at  P^ 

DENOM  -  the  denominator  of  equations  (78)  ^nd  (8l) 

XB  -  X  coordinate  of  P^ 

YB  -  Y  coordinate  of  P 

DENO  -  the  denominator  of  equations  (78)  and  (8l)  used  to  calculate 
XB  and  YB 

XL  -  X  coordinate  of  P2 
YL  -  Y  coordinate  of  P 

2 

DEN  -  the  denominator  of  equations  (78)  and  (8l)  used  to  calculate 
XL  and  YL 
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DXB  -  the  absolute  difference  in  the  X  coordinates  of  P„  and  P, 

0  1 

DYB  -  the  absolute  difference  in  the  Y  coordinates  of  P  and  P 

0  1 

DYB 

AB  •  the  orctan  of 

DXL  -  the  absolute  difference  in  the  X  coordinates  of  P^  and  P„ 

0  2 

DYL  -  the  absolute  difference  in  the  Y  coordinates  of  P^  and  P^ 

DYL 

AL  -  the  arc  tan  of  =rrs- 


DMM  -  the  distance  on  the  projection  plane  between  PQ  and  P^ 

DPM  -  the  distance  on  the  projection  plane  between  PQ  and  Po 

DME  -  distance  on  the  datum  surface  between  P  and  P 

0  1 

DPE  -  the  distance  on  the  datum  surface  between  P  and  P 

0  2 

The  following  is  a  brief  explanation  of  the  Source  Language 
Statements  of  the  program  included  at  the  end  of  this  appendix: 

1  -  35  J  set  up  the  projection  parameters  and  are  self  explanatory 
when  considered  with  the  term  ceTinitions  given  above 
36  -  45  :  calculate  the  X  and  Y  coordinates  of  PQ 
46  -  71  :  check  the  location  of  PQ  to  determine  if  it  is  visible 
from  the  projection  center 
72  -  80  :  calculate  the  X  and  Y  coordinates  of  P^ 

8l  -  87  J  calculate  the  X  and  Y  coordinates  of  P 

2 

88  -  105  t  calculate  the  angular  distortion  in  the  projected 
intersection  of  a  parallel  and  a  meridian 
106  -  ll6  :  calculate  the  meridian  and  parallel  linear  distortions 
117  -  153  •  dictate  the  size  and  shape  of  grid  to  be  computed 

154  -  159  J  calculate  the  longitude  values  that  will  be  used  for 

column  headings  in  the  grid 

9L 


160  •  167  *  print  out  the  projection  parameters 
l68  -  197  I  print  eut  the  numerical  grid 

The  program  contained  on  the  next  seven  pages  vas  used  t 
produce  Example  5  of  this  thesis. 
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